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5^ : 1. INTRODUCTION. 

These notes are an expanded version of the lectures I gave at the 1991 ICTP 
Spring School on String Theory and Quantum Gravity. Here I have attempted to 
>• ■ review, from my own personal viewpoint, some of the exciting developments in two- 
dimensional string theory that have taken place over the last year and a half. Because 
of the multitude of new results, and since the field is still developing rapidly, a com- 
prehensive review must await a later date. These notes are mainly devoted to the 
matrix model approach 1 that has truly revolutionized the two-dimensional Euclidean 
quantum gravity. Recently this approach has led to the exact solution of quantum 
gravity coupled to conformal matter systems with c < 1. 2 ~ 6 

These lectures are about the c = 1 model 5,6 that is both the richest physically 
and among the most easily soluble. This model is defined by the sum over geometries 



★ Lectures delivered at the ICTP Spring School on String Theory and Quantum Gravity, Trieste, 
April 1991. 
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embedded in 1 dimension. The resulting string theory is, however, 2-dimensional 
because the dynamical conformal factor of the world sheet geometry acts as an extra 
hidden string coordinate. 7 In fact, this is the maximal bosonic string theory that 
is well-defined perturbatively. If the dimensionality is increased any further, then a 
tachyon appears in the string spectrum and renders the theory unstable. The matrix 
model maps the theory of surfaces embedded in 1 dimension into a non-relativistic 
quantum mechanics of free fermions, 8 from which virtually any quantity can be 
calculated to all orders in the genus expansion. The second-quantized field theory 
of non-relativistic fermions can be regarded here as the exact string field theory. A 
transformation of this theory to a rather simple interacting boson representation 9 ' 10 
will also be given. I will argue that the 2-dimensional string theory is the kind 
of toy model which possesses a remarkably simple structure, and at the same time 
incorporates some of the physics of string theories embedded in higher dimensions. 
The simplicity of the theory is apparent in the matrix model approach, but much of 
it remains obscure and mysterious from the point of view of the continuum Polyakov 
path integral. Until we develop a better insight into the "miracles" of the matrix 
models, we may be sure that our understanding of string physics is very incomplete. 

These notes are mainly a review of published papers, but I also included a few 
new observations and results. In section 2 I review the formulation of quantum grav- 
ity coupled to a scalar field as a sum over discretized random surfaces embedded in 1 
dimension. I will show that the matrix quantum mechanics generates the necessary 
statistical sum. In section 3 I exhibit the reduction to free fermions and define the 
double-scaling limit. 5 The sum over continuous surfaces is calculated to all orders in 
the genus expansion. In section 4 I show how to calculate the exact string amplitudes 
using the free fermion formalism. Section 5 is devoted to the continuum path integral 
formalism 11 where the Liouville field acts as an extra dimension of string theory. 
Some of the exact matrix model results will be reproduced, but this approach still 
falls far short of the power of the matrix model. In section 6 I discuss the special 
states which exist in the spectrum in addition to the massless "tachyon". These 
states, occurring only at integer momenta, are left-overs of the transverse excitations 
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of string theory. Remarkably, they generate a chiral W^i+oo algebra. In section 7 I 
use the matrix model to formulate the exact string field theory both in the fermionic 
and in bosonic terms. I present some manifestly finite bosonic calculations which are 
in perfect agreement with the fermionic ones. In section 8 I discuss the new physical 
effects that arise when the random surfaces are embedded in a circle of radius R: the 
it! — > 1/R duality and its breaking due to Kosterlitz-Thouless vortices. 6 This model 
can be regarded either as a compactified Euclidean string theory or as a Minkowski 
signature string theory at finite temperature. In section 9 I show that in the matrix 
quantum the effects of vortices are implemented by the states in the non-trivial repre- 
sentations of SU (iV). 6 ' 12 Including only the SU (N) singlet sector gives the vortex-free 
continuum limit where the sum over surfaces respects R — > 1/R duality. 6 In section 
10 I use the thermal field theory of non-relativistic fermions to find exact amplitudes 
of the compactified string theory. 13 Finally, in section 111 use the matrix chain model 
to solve a string theory with a discretized embedding dimension. 6 Remarkably, when 
the lattice spacing is not too large, this theory is exactly equivalent to string theory 
with a continuous embedding. 



In this section I will introduce the discretized approach to summation over random 
surfaces embedded in one dimension and its matrix model implementation. If we 
parametrize continuous surfaces by coordinates (a\, (72), then the Euclidean geometry 
is described by the metric g^ v (a), and the embedding - by the scalar field X(a). 
Thus, the theory of random surfaces in one embedding dimension is equivalent to 2-d 
quantum gravity coupled to a scalar field. In the Euclidean path integral approach to 
such a theory, we have to sum over all the compact connected geometries and their 
embeddings, 



The integration measure is defined modulo reparametrizations, i.e. different func- 



2. DISCRETIZED RANDOM SURFACES 
AND MATRIX QUANTUM MECHANICS. 




topologies 



(2.1) 
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tions describing the same geometry-embedding in different coordinates should not be 
counted separately. We assume the simplest generally covariant massless action, 




d 2 °V9( —^dpXdvX + $R + 4A) . 



(2.2) 



In 2 dimensions the Einstein term is well-known to give the Euler characteristic, the 
topological invariant which depends only on the genus of the surface h (the number 
of handles) , 



genus h surfaces can be thought of as a diagram of string theory with h loops, we 
identify the string coupling constant as go = e*. 

The main problem faced in the calculation of the Euclidean path integral of eq. 
(2.1) is the generally covariant definition of the measure for the sum over metrics 
[Dg^v] and of the cut-off. One may attempt to do this directly in the continuum. 
Considerable success along this route has been achieved when metrics are described in 
the light-cone gauge, 14 or in the more traditional conformal gauge. 15 We will carry 
out some comparisons with the continuum approach in the course of the lectures. 

The main subject of this section is a different approach to summing over geome- 
tries, which has so far proven to be far more powerful than the continuum methods. 
In this approach one sums over discrete approximations to smooth surfaces, and then 
defines the continuum sum by taking the lattice spacing to zero. 1 We may, for in- 
stance, choose to approximate surfaces by collections of equilateral triangles of side a. 
A small section of such a triangulated surface is shown in fig. 1. The dotted lines show 
the lattice of coordination number 3, which is dual to the triangular lattice. Each face 
of the triangular lattice is thought of as flat, and the curvature is entirely concentrated 
at the vertices. Indeed, each vertex I of the triangular lattice has a conical singularity 
with deficit angle f (6 — qj), where qi is the number of triangles that meet at I. Thus, 




(2.3) 



Thus, the weighting factor for a surface of genus h is (exp <fr) 2h 2 . Since the sum over 
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Fig. 1. A small section of triangulated surface. Solid lines denote the trian- 
gular lattice A, and dotted lines - the dual lattice A. 

the vertex I has a 5-function of curvature with positive, zero, or negative strength 
if qi < 6, =6, > 6 respectively. Such a distribution of curvature may seem like a 
"poor man's version" of geometry. However, in the continuum limit the size of each 
face becomes infinitesimal, and we may define smoothed out curvature by averaging 
over many triangles. In this way, the continuum field for the metric should appear 
similarly to how the quantum field description emerges in the continuum limit of the 
more familiar statistical systems, such as the Ising model, the XY model, etc. Later 
on, we will show some strong indications that the discretized approach to summing 
over geometries is indeed equivalent to the continuum field-theoretic approach. 

The essential assumption of the discretized approach is that sum over genus h 
geometries, jlDg^]^, may be defined as the sum over all distinct lattices A, with the 
lattice spacing subsequently taken to zero. f[DX] is then defined as the integral over 
all possible embeddings of the lattice A in the real line. For simplicity, the lattices 
A may be taken to be triangular, but admixture of higher polygons should not affect 
the continuum limit. In order to complete the definition of the discretized approach, 
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we need to specify the weight attached to each configuration. This can be defined by 
simply discretizing the action of eq. (2.2) and counting each distinct configuration 
with weight e~ s . We will find it convenient to specify the embedding coordinates 
X at the centers of the triangles, i.e. at the vertices i of the dual lattice. Then the 
discretized version of J (Pa^fgg^d^XdpX is ~ J2<ij>(Xi — Xj) 2 , where the 
sum runs over all the links < ij > of the dual lattice. Similarly, 

/ #oy/g - J«V (2.4) 

where V is the number of triangles, or, equivalently, the number of vertices of the 
dual lattice A. Thus, the discretized version of the path integral (2.1) is 

V r 

2(go, «) = £ E * V II / dX * II G ^ - X i) > ( 2 - 5 ) 

h A i=l (ij) 

where A are all triangular lattices of genus h, k ~ exp(— \/3Xa 2 /Atc), and (for some 
choice of a') G(X) = exp(— \X 2 ). If the continuum limit of eq. (2.5) indeed describes 
quantum gravity coupled to a scalar field, then there should exist a whole universality 
class of link factors G, of which the gaussian is only one representative, that result in 
the same continuum theory. 

A direct evaluation of the lattice sum (2.5) seems to be a daunting task, which 
is still outside the numerical power of modern computers. Fortunately, there exists 
a remarkable trick which allows us to exactly evaluate sums over surfaces of any 
topology using only analytic tools. As was first noted by Kazakov and Migdal, 16 a 
statistical sum of the form (2.5) is generated in the Feynman graph expansion of 
the quantum mechanics of a N x N hermitian matrix* Consider the Euclidean path 



★ Similar tools work for other simple matter systems coupled to two-dimensional gravity. For 
example, as discussed in other lectures in this volume, in the case of pure gravity the sum over 
discretized surfaces is generated simply by an integral over an N x N hermitian matrix. 
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integral 



D N ~$(x)exp 




(2.6) 



where x is the Euclidean time and U = 2^7$ 2 — ^t^ 3 - The parameter (3 enters as the 



3! 



inverse Planck constant. By a rescaling of $ eq. (2.6) can be brought to the form 

T/2 



/ 



D^'^aOexp 




2a' 3! 



(2.7) 



where k = yjN/(3 is the cubic coupling constant. The connection with the statistical 
sum (2.5) follows when we develop the graph expansion of Z in powers of k. The 
Feynman graphs all have coordination number 3, and are in one-to-one correspon- 
dence with the dual lattices A of the discretized random surfaces (fig. 1). The lattices 
A dual to the Feynman graphs can be thought of as the basic triangulations. One 
easily obtains the sum over all connected graphs In Z, 



h A i=lj , 00 (ij) 



-\xi-Xj\/a' 



(2.8) 



This is precisely of the same form as eq. (2.5) which arises in two-dimensional quan- 
tum gravity! The Euclidean time x assumes the role of the embedding coordinate 
X in eq. (2.5). We note that the role of the link factor G is played by the one- 
dimensional Euclidean propagator. Only for this exponential G can we establish the 
exact equivalence with the matrix model. This does not pose a problem, however, as 
we will find plenty of evidence that the continuum limit of the model (2.5) with the 
exponential G indeed describes quantum gravity coupled to a scalar, eq. (2.1). It is 
evident from (2.8) that the parameter a' sets the scale of the embedding coordinate. 
In fact, we have normalized a' so that in the continuum limit it will precisely coincide 
with the definition in eq. (2.2). Whenever a' is not explicitly mentioned, its value 
has been set to 1. 
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Further comparing eqs. (2.8) and (2.5), we note that the size of the matrix 
N enters as 1/go- Let us show why. Each vertex contributes a factor ~ N, each 
edge (propagator) ~ 1/N, and each face contributes N because there are as many 
index loops as there are faces. Thus, each graph is weig hted by N V ~ E+F which, 
by Euler's theorem, equals jV 2 ~ 2ft . The expansion of the free energy of the matrix 
quantum mechanics in powers of 1/N 2 automatically classifies surfaces according to 
topology. It would seem that, in order to define a theory with a finite string coupling, 
it is necessary to consider finite N which, as we will see, is associated with severe 
difficulties. Fortunately, this naive expectation is false: in the continuum limit the 
"bare" string coupling 1/N becomes infinitely multiplicatively renormalized. Thus, if 
N is taken to oo simultaneously with the world sheet continuum limit, then we may 
obtain a theory with a finite string coupling. This remarkable phenomenon, known as 
the "double-scaling limit" 2 ' 5 , will allow us to calculate sums over continuum surfaces 
of any topology. 

In order to reach the continuum limit, it is necessary to increase the cubic coupling 
k = \/N/(3 ~ N° until it reaches the critical value k c where the average number of 
triangles in a surface begins to diverge. In this limit we may think of each triangle as 
being of infinitesimal area \/3a 2 /A, so that the whole surface has a finite area and is 
continuous. In the continuum limit k — > k c we define the cosmological constant 

A = 7r(/t 2 - k 2 ) . (2.9) 

Recalling that k = K, c exp(—\^3\a 2 /An), we establish the connection between A and 
the physical cosmological constant A: A ~ A/a 2 . 

Above we have sketched the connection between matrix quantum mechanics and 
triangulated random surfaces. It is easy to generalize this to the case where, in 
addition to triangles, the surfaces consist of other n-gons: we simply have to add to 
the matrix potential U (<&) monomials $ n , n > 3. For consistency, the continuum limit 
should not be sensitive to the precise manner in which the surfaces are discretized. 
In the next section we will give the exact solution of the matrix model and show that 
the continuum limit is indeed universal. 
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3. MATRIX QUANTUM MECHANICS AND FREE FERMIONS. 



In the previous section we established the equivalence of the sum over connected 
discretized surfaces to the logarithm of the path integral of the matrix quantum 
mechanics. On the other hand, in the hamiltonian language, 

Z =< f\e~ /3HT \i> . (3.1) 

Thus, as long as the initial and final states have some overlap with the ground state, 
the ground state energy Eq will dominate the T — > oo limit of eq. (3.1), 

lim ^ = _ PEq . (3.2) 

The divergence of the sum over surfaces proportional to the length of the embedding 
dimension arises due to the translation invariance. In order to calculate the sum 
over surfaces embedded in the infinite real line R 1 , all we need to find is the ground 
state energy of the matrix quantum mechanics. To this end we carry out canonical 
quantization of the SU (N) symmetric hermitian matrix quantum mechanics. 8 ' 17 The 
Minkowski time lagrangian 

L = Tr{±$ 2 - [/($)} (3.3) 

is symmetric under time-independent SU(N) rotations $(t) — > V^(t)V. It is con- 
venient to decompose $ into N eigenvalues and N 2 — N angular degrees of freedom, 

= &{t)A{t)Q{t) (3.4) 

where f2 e SU(N) and A = diag{\\, A2, . . . , Ajv)- The only term in the lagrangian 
which gives rise to dependence on Q is 

Tr <i> 2 = Tr A 2 + Tr[A, QQ^] 2 . (3.5) 

To identify the canonical angular degrees of freedom, the anti-hermitian traceless 
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matrix can be decomposed in terms of the generators of SU(N), 



N-l 

X ^ aiHi (3.6) 



V i<j i—1 



l <3 

where Hi are the diagonal generators of the Cartan subalgebra, and the other gener- 
ators are defined as follows. T y - is the matrix M such that M y - = Mji = 1, and all 
other entries are 0; T y - is the matrix M such that My = —Mji = —i, and all other 
entries are 0. Substituting eq. (3.6) into eqs. (3.5) and (3.3), we find 

L = J2 + u &)) + i - A i) 2 (4 + ■ ( 3 - 7 ) 

i i<j 

In deriving the hamiltonian it is crucial 8 that the measure of integration expressed 
in terms of £1 and A assumes the form £>$ = VQ^dXiA 2 ^), where A(A) is the 
Vandermonde determinant rL<j(^i ~~ ^j)- Because of the non-trivial Jacobian, the 
kinetic term for the eigenvalues becomes 

1=1 

It is left as an exercise for the reader to show that this can be rewritten as 
It follows that the total hamiltonian is 



» - "TO) E + £ UM + g g±J , (3,0) 

where Hy and Ily are the momenta conjugate to cty and respectively, i.e., they are 
the generators of left rotations on Q, Q — > Furthermore, since L is independent 
of aii, the wave functions must obey the constraints = 0. These constraints 

arise because transformations VI — > AQ, where A is a diagonal SU(N) matrix, do not 
change <£>. 
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The constraints require that only those irreducible representations of SU(N) 
that have a state with all weights equal to zero are allowed in the matrix quan- 
tum mechanics. 12 For instance, the fundamental representations are excluded, and 
the simplest non-trivial representation is the adjoint. It is not hard to classify all the 
irreducible representations allowed in the matrix quantum mechanics, i.e., the ones 
that have states with zero weights. In the standard notation, the Young tableaux 
of such representations consist of i + j columns, with the numbers of boxes in the 
columns (from left to right) 



where n\ < 112 < . . . < rii < N/2 and mi < . . . < m j < N/2. 

Since the angular kinetic terms in (3.10) are positive definite, it is clear that we 
have to look for the ground state among the wave functions that are annihilated 
by Uij and ITy, i.e. in the trivial (singlet) representation of SU(N). The singlet 
wave functions are independent of the angles Q, and are symmetric functions of the 
eigenvalues, Xsym(A). Because of the special form of the hamiltonian (3.10), the 
eigenvalue problem Hxsym = Exsym assumes the remarkable form 8 



where ^(A) = A(\)x S ym(^) is by construction antisymmetric. Since the operator 
acting on \1> is simply the sum of single-particle non-relativistic hamiltonians, the 
problem has been reduced to the physics of iV free non-relativistic fermions moving 
in the potential U(\). 8 This equivalence to free fermions is at the heart of the exact 
solubility of string theory in two dimensions. 

The iV-fermion ground state is obtained by filling the lowest N energy levels of 



(N - mi, N - m 2 , . . . , N - mj,rii, . . . ,n 2 , m) 



(3.11) 




(3.12) 



hi = 



2f3 2 d\i 



2 + U(\i) , 
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h, 

N 

E = J2 £ i- ( 3 - 13 ) 

i=l 

The potential arising from triangulated random surfaces is U (A) = 5A 2 — -^jA 3 , which is 
unbounded from below and does not seem support any stable states. Recall, however, 
that in order to find the genus expansion of the sum over surfaces, we are only 
interested in the expansion about the classical limit f3 — > 00 (h — > 0) in powers of 
1/ ' 1 . Classically, there is a continuum of bounded orbits to the left of the potential 
barrier. Semiclassically, these orbits are associated with energy levels (fig. 1) whose 
spacing is 0(1/ (3), and whose decay time is exponentially long as (3 — > 00. This 
instability cannot be seen in the context of the semiclassical expansion in powers of 
1/ ' 1 . Our goal, therefore, is to find the expansion of the ground state energy (3.13), 
where are the semiclassical energy levels. 

To leading order, e n are determined by the Bohr-Sommerfeld quantization condi- 
tion, 

/2tt 
p n (\)d\ = — n (3.14) 

where p n (A) = \/2(e n — U (A)), and the integral is over a closed classical orbit. In 
particular, the position of the Fermi level fip = ejy is determined by 

A+ 

J y/2fr F -U(\))=1Tj, (3.15) 

where A_ and A + are the classical turning points. Eq. (3.15) shows that jip is an 
increasing function of N//3. Clearly, both fip and Eq become singular at the point 
where N/(3 — > such that jip equals the height of the barrier fx c . Recall that the 
cubic coupling of the matrix quantum mechanics is k — a/ N/(3, and that we expect 
to find a singularity in the sum over graphs when k becomes large enough that the 
average number of vertices in a Feynman graph (or equivalently, of triangles in a 
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Fig. 2a) N fermions in the Fig. 2b) The double-scaling 

asymmetric potential arising limit magnifies the quadratic 

directly from the triangulated local maximum, 
surfaces. 



random surface) begins to diverge. In the vicinity of this singularity the continuum 
limit of quantum gravity can be defined. Now we have identified the physics of this 
singularity in the equivalent free fermion system: it is associated with spilling of 
fermions over the barrier. 

We are interested in the singularity of Eq as a function of the cosmological con- 
stant A = ir{K 2 c — ^). The parts of E(A) analytic in A are not universal and can be 
dropped. Indeed, if we calculate the sum over surfaces as a function of the fixed area 
by inversely Laplace transforming E(A), the analytic terms give contributions only 
for zero area. 

Let us introduce the density of eigenvalues 

" n 

in terms of which 
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liu^^A = (3E = (3 2 J p(e)ede 



(3.17) 



o 



Differentiating the first equation, we find 

<9A 

— = np(fj, F ) , (3.18) 

which can be integrated and inverted to obtain fx = fi c — fip as a function of A. 
Differentiating the second equation, we find |^ = ^f3(fx — fx c ). After addition to E 
of an irrelevant analytic term, E — > E + ^(3fx c A, we find 

m = -^ (3 - 19 > 

These equations show that the continuum limit of the sum over surfaces is determined 
by the singularity of the single particle density of states near the top of the barrier. 
Using the WKB approximation, it was found in ref. 16 that 

A+ 

If d\ 1 , 2 



P^f) = ~ j= ~--]nn + 0(l//3*) (3.20) 



It follows that A = -fx\a.fx + <D(l/f3 2 ), and 

The critical behavior of p comes from the part of the integral near the quadratic 
maximum of the potential. In the continuum limit fx — > 0, the classical trajectory at 
[xf spends logarithmically diverging amount of time near the maximum. Therefore, 
the WKB wave functions are peaked there. In fact, we will show that, at any order 
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Fig. 3a) N fermions in a Fig. 3b) The system rel- 

double well potential. evant to the double-scaling 

limit. 



in the WKB expansion, the leading singularity of p is entirely determined by the 
quadratic nature of the maximum. To see this, let us rescale the coordinates \/]3(\ — 
A c ) = y, and e = /i c — e. The single particle eigenvalue problem becomes 5 

+ \y 2 + 0(i/V3)v 3 ) m = P*Kv) (3-22) 

For a finite f3e, all the terms beyond the quadratic are irrelevant in the limit (3 — > oo 
because they are suppressed by powers of (3. This demonstrates the universality of 
the continuum limit: any potential possessing a quadratic maximum will yield the 
same sum over continuous surfaces. In effect, our rescaling has swept all the non- 
universal details of the potential out to infinity. This implies, for instance, that, 
instead of working with the potential of fig. 2, we can use the double well potential 
U(X) = — \\ 2 + A 4 (fig. 3). After rescaling y/~j3\ = y, the Schrodinger equation again 
reduces to the motion in the inverted harmonic oscillator. The only difference is that 
now fermions fill both sides of the maximum. To all orders of the WKB expansion 
this simply multiplies the density of states, and therefore the free energy, by a factor 
of 2. The double well potential does not suffer from the non-perturbative instabilities: 
the ground state energy can in principle be calculated even for finite N and (3. We 
will often adopt the symmetric case because it is more convenient for calculations, 
remembering that we have to divide by 2 to find the sum over triangulated surfaces. 
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From eq. (3.21) we know the sum over all surfaces of spherical topology embedded 
in R 1 . In order to find the sums over surfaces of any topology, we need to develop a 
systematic WKB expansion of the ground state energy. According to eqs. (3.18) and 
(3.19), all we need to know is the complete expansion of the single particle density of 
states p(/i). If we consider /i such that (3fi is Oil), then this problem has a remarkably 
simple solution, 5 because eq. (3.22) shows that in this energy range the system is 
purely quadratic. Let us write the density of states as 



p(li) =Tr 5(h + fa) = -ImTr 

7T 

1 d 2 1 2 

h0 = -2d^-2 y ■ 



ho + fifj, — it 



(3.23) 



We will evaluate the trace in position space in order to take advantage of the fact 
that the resolvent of the hamiltonian of a simple harmonic oscillator is well known, 

oo Vf 

(V,\ n, ;4l ,, r\Vi)= I iT^ T I Vy(t)e- I>i<M'^ 



= [ dTe-^T / u c -Lo[(y 2 f +y?)coshLoT-2y m ]/2 S mhcuT 

J V 27rsinhc;T 



In our case, the frequency uj is actually imaginary (we have an inverted harmonic 
oscillator), so we must define the resolvent by analytic continuation. We rotate uj — > 
—i while simultaneously rotating T — > iT, and the relevant resolvent is 



■oo 



/ i \ 1 )=l \ dTe -W 4(^ 2 )coshr-2^]/2si n hr 

yyn h + P/i-ie m/ J V 2vrsinhT 



(3.25) 

Now it is easy to derive the derivative of the density of states 5 

oo oo 

dp 1 9 f J ,, 1 i \ 1 T _ f Jrr , —iQuT T/2 



-Im / dy(y\-r—± -\y) = -Im / dTe~^ T - 



d(Pn) Trd(Pn) J *^ l h + Ppi-ie*"' vr J sinh(T/2) 

-oo 

(3.26) 

Note that only for the purposes of deriving the large expansion can we reduce 
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the problem to pure inverted harmonic oscillator. The integral representation (3.26) 
is not designed to correctly give the non-perturbative terms (9(e _/3 ^). In fact, such 
terms depend on the details of the potential away from the quadratic maximum, and 
are not universal. While each term in the large expansion is related to the sum 
over all geometries of a certain genus, we do not yet fully understand the geometrical 
or physical meaning of the non-perturbative terms. 

After expanding eq. (3.26) and integrating once, we find the complete series of 
corrections to the leading order estimate in eq. (3.20), 

p(^) = -(-ln/i+f;(2 2 - 1 -l)^(2^r 2 -l . (3.27) 

71 { m=l m J 

In solving eq. (3.18) for /i in terms of A, it is useful to introduce parameter /iq 
defined by 

A = -/ioln/io = /io| hi/!o| • (3.28) 

Then, 

Ax = /xo{l-^-f; (2 2 - 1 ~ 1) 'f 2m ' - ; (2^o)" 2m + O (*) ) . (3.29) 
[ ln ^°m^l m(2m-l) Vhi>o/J 

Integrating (3.19) we finally arrive at the complete genus expansion of the sum over 
surfaces 

-PE = ±- ((2^ ) 2 ln / ,o - \ ln^o + £ ^tSt^ ^^] " 
87T 1 3 m(m + l)(2m + 1) J 

(3.30) 

Clearly, the sensible way to approach the continuum limit is through "double scaling" 
(3 — > oo and /iq — > in such a way that = l/g s t is kept fixed. 5 ' 6 Thus, quite 
remarkably, we have found closed form expressions for sums over continuous surfaces 
of any topology embedded in one dimension. The sum over genus h surfaces, given 
by the coefficient of g^~ 2 , grows as (2h)\. This behavior is characteristic of closed 
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string theories, while in field theory we typically find that the sum over /i-loop 
diagrams exhibits a slower growth ~ h\. The badly divergent perturbation series 
(3.30) indicates that the theory is not well-defined non-perturbatively. 

Eq. (3.30) shows explicitly how the bare string coupling 1/(3 gets multiplicatively 
renormalized. Compared to the c < 1 matter coupled to gravity, 2-4 the new feature 
here is that the renormalization is not simply by a power of the cosmological constant 
A. It appears that A itself is multiplicatively renormalized: A — > /j,q A/|lnA|. 
Another peculiarity of the c = 1 matter is that the sums over spherical and toroidal 
surfaces diverge in the continuum limit as | In /xo | - In section 5 we will argue that all 
these features have a natural explanation in the context of Liouville theory. Before 
we do that, however, we further exhibit the power of the matrix model by calculating 
some correlation functions to all orders in the topological expansion. 

4. CORRELATION FUNCTIONS. 

The simplest kind of matrix model operator is Tr$ n (:r) where n is any finite 
integer. This operator inserts a vertex of order n in all possible places on the dual 
graph A. In terms of the basic lattice A, this amounts to cutting an n-gon hole in 
the surface, pinning it at the embedding coordinate x, and integrating the position 
of the hole all over the surface. In the continuum limit, for any finite n the size of 
the hole becomes infinitesimal, and it is commonly referred to as "a puncture". The 
continuum expression for the operator above should be roughly 

J d 2 a^5{X(a)-x) . (4.1) 

By Fourier transforming, we derive the continuum translation of a related operator 

J dxe iqx Tr$ n (x) -> J d 2 a^je iqX , (4.2) 

which resembles the basic "tachyon" operator of string theory. Below we outline the 
procedure for calculating the correlation functions of such operators in the matrix 
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model. In section 5 we show that the results can indeed be reproduced in the string 
theoretic formalism. 

Since the matrix model has been reduced to a system of free fermions, it is con- 
venient to perform calculations in the formalism of non-relativistic second quantized 
field theory. The second quantized fermion field is defined as 

oo 

$(\,x)= J due iux a e (u)^ e (u,\) , (4.3) 

— oo 

where ip e (v, A) are the single fermion wave functions of energy v, and the associated 
annihilation operators a e {v) satisfy {aj(i/), a e (v')} = 5 e ^'5{i> — v'). Setting A = \ply 
in eq. (3.22), we find that ip e (u, A) are solutions of the Whittaker equation 

+ *)1> = ^- (4-4) 

and can therefore be expressed in terms of Whittaker functions W(u, A). 19 Their 
spectrum is continuous and doubly degenerate. Here e denotes the parity of the 
wavefunction ip e (v, A) and repeated e are summed over. 

The Euclidean second-quantized field theory with chemical potential fifi is defined 
by the action 

oo oo 2 2 

S= I dxj dxtf(-± + ±- + ^ + Mt (4.5) 



-oo — oo 



where x is the Euclidean time. The ground state of the system (fig. 3b) satisfies 

a e (v)\(3fi) = 0, v<Ppt, 

(4.6) 

al(v)\P/j) = 0, v > fifx . 

Consider correlation functions of operators of the type 

0(q) = J dxe iqx Trf($(x)) , (4.7) 
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where / is any function. In the second quantized formalism, they translate into 



0(q) = J dxe iqx J d\f(\)fl$(x,\) . 



(4.8) 



Thus, to find the connected correlation functions of any set of such operators, we 
need to calculate 

n „ 

G( gi ,Ai;...;g n ,A n )=JI / dx i e i * Xi (P»\fl$(x 1 ,\ 1 )..4U(x n ,\ n )\P f i) c . (4.9) 

i=i 

Application of Wick's theorem reduces this to the sum over one-loop diagrams with 
fermion bilinear insertions in all possible orders around the loop. A convenient formula 
for the Euclidean Green function is 19 

oo 

S E (x h Ai; x 2 , A 2 ) = J due-^ l3 ^ Ax {9(Ax)9(u - - 0(-Ax)9(Pfx - is)}x 



— oo 
oo 



oo sgn(p)oo 
dp -ipAx 
' ' 27T 6 



p + i(v — (3fi) 

)oo 

J dse~ ap+i ^ s {\i\e 2isH \\ 2 ) 







(4.10) 



where Ax — x\ — x 2 , and we have used 19 



oo 

/ 



dve- ws ^{yMW(y, A 2 ) = (Ai|e 2iaff |A 2 > 



= exp — 
V — 47ri sh s V 4 



A 1 + A 2 2AiA 2 
th s sh s 



(4.11) 



with if = l/2p 2 — 1/8A 2 . With the Green function (4.10), after integration over the 
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loop momentum, Moore derived the representation 



oo eioo e„_ioo 

-^G(q 1 ,\ 1 ;...;q n ,\ n )=t n+1 5(J2u) E J ^ J d Sl . . . J ds n _ 



(4.12) 

where Q£ = + ... + q a (k)> an d e fc — s 9 n [Qk\- 111 principle, any correlation 

function of operators of type (4.7) can be found by integrating G. A convenient way 
of calculating correlations of puncture operators is to first introduce the operator that 
creates a finite boundary of length /, with momentum q injected into the boundary 19 



0{l,q)= J dxe iqx Tre- l *W 



(4.13) 



The puncture operator should be the leading term in the small / expansion of 0(1, q). 
The correlations of puncture operators P(q) can be extracted according to 



(4.14) 



The details of the calculations are highly technical and can be found in ref. 19. Here 
we will simply state the results for the two, three, and four point functions. 



(P(q)P(k)) = -5(q + k)[T(l-\q\)] 2 ^ 



1 



1 



(4.15) 



(P(qi)P(q2)P(qs)) = S(J2 *) II ( r (* ~ l«*l)A* 

i=l 

4 



[M\ 1 



+ nCksj - r) U(q\ + q$) + lOqfql - 10(g? + g 2 2 )|g 3 | - 5(qi ~ q2? + 12|g 3 | + l) 



X 



5760(/fy) : 



+ . . . 



(4.16) 
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where qi, q2 > 0. 



(P(9i)P(ffi)P(ffl)P(?4)> = fa 1 " 1^1 V ¥ ) x 

i=i (4.17) 

[(191 + © I + 1 91 + ©| + |gi + 94| - 2) + 0(l/(/fy) 2 )] . 

The new feature of the c = 1 correlators, compared to theories with c < 1, is that they 
have genuine divergences occurring at quantized values of the external momenta. 20 In 
this sense, the c = 1 model is the most similar to critical string theory where the 
amplitudes are well-known to have divergences associated with the production of on- 
shell physical particles. The remarkable feature of the c = 1 amplitudes is that all 
the divergences factorize into external leg factors. 

Because of the divergences, eqs. (4.15)-(4.17) are strictly valid only if none of 
qi are integer. If any of the momenta have integer values, then additional terms 
need to be taken into account, which leads to cancellation of the infinity. 20 The 
physical meaning of this is simple: in the matrix model there is an explicit ultravi- 
olet and infrared cut-off, and therefore there are no genuine zeroes or poles of the 
correlation functions. Instead, they are regularized according to 1/0 — > |ln/xo|, and 
— > l/|ln//o|- Note that in eqs. (4.15)-(4.17) fx can be replaced by liq. Thus, the 
correlation functions scale as powers of /iq, and not as powers of A found in all the 
c < 1 theories coupled to gravity. 2-4 This provides further evidence that, for c = 1, 
a renormalization of the cosmological constant takes place. In the next section we 
sketch a Liouville theory explanation of some of the new effects peculiar to c = 1. 
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5. LIOUVILLE GRAVITY COUPLED TO c = 1 MATTER. 



In this section we outline the continuum approach to evaluating the sum over 
surfaces (2.1). In order to integrate over all metrics, in 2 dimensions we may pick the 
conformal gauge = g^vij)^^ where g^ v depends on a finite number of modular 
parameters collectively denoted as r. For spherical topology, there are no moduli, 
and we may choose g^ u = 8^ v . Eq. (2.2) shows that classically the Liouville field <fi is 
non-dynamical. It is well-known, however, 11 that the kinetic term for is induced 
through quantum effects arising from dependence on <p of the path integral measure. 
Additional modification of the dynamics of <fi has been proposed by David and by 
Distler and Kawai, 15 who argued that we can simplify the integration measure at the 
expense of further renormalization of the action for such that the ghost, X and (f> 
fields combine into a conformally invariant theory with net conformal anomaly equal 
to zero. Following their approach, we assume that the original path integral (2.1) is 
transformed in the conformal gauge into 

Z = / [dr^DbDc^D^DX^- 3 , 

* (5.1) 
S = — / tfoyfl (d^Xd^X + d^(j) - 4R<f) + 4Ae~ 2< ^ , 

and Si iC is the standard free action for the b and c ghosts. The exponential interaction 
term originates from the cosmological term in the action (2.2).* 

With the action (5.1), the simplest set of conformally invariant operators are 

T(q) = J- / d 2 <j^fg~e iqX+ <^ (5.2) 

7TP J 

These operators inject embedding momentum q into the world sheet. They should be 
thought of as the Liouville theory implementation of the operators (4.2). In a theory 

* Polchinski has argued 21 that, for c = 1, the interaction term is not simply an exponential. 
We will instead work with the exponential interaction, and will show that it is also possible 
to explain correlation functions within this framework. Our argument docs not substantially 
differ from Polchinski's, but we prefer working with the exponential because it was shown by 
Curtwright and Thorn 22 that with this form of interaction the conformal invariance can be 
maintained exactly in the continuum limit. 
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with A = 0, T(q) is conformally invariant for e(q) = —2 ± \q\. For A > 0, however, 
it was shown in ref. 23, 24 that only the operator with e(q) = — 2 + |g| exists. This 
conclusion agrees with the scaling of correlation functions in the matrix model. Let 
us now use the Liouville formalism to calculate the spherical correlation functions. 

As discovered in ref. 25, in performing the path integral it is convenient to 
decompose <p = <f)Q + <p, and integrate first over the zero mode 0o- Then, since the 
interaction is purely exponential, we obtain 



ri T ^)) =#- 2 '£>)i(£) r(- a ) 
1=1 / h ' 

J [dr][DX][D^][DbDc]e- s ^- So (^J d 2 a^ e - 2 ^ U^p J ^ a ^ e 



s 

~Jq t X+(-2+\ qi \)4> 



(5.3) 



where 



s ° = h J d2ay ^ ( d ^ d ^ + d ^i) ( 5 - 4 ) 

and s = 2(1 — h) — N + \ Y^i=\ 1^1- The factor r(— s) indicates that correlation 
functions with non- negative integer s have a divergence from the volume of 0o- This 
explains the divergences in the 0-point functions on the sphere and torus found in 
eq. (3.30). Now, the difficulty is to deal with the insertion of T(0) s . Goulian and 
Lee 26 proposed to calculate eq. (5.3) first for integer s, where the rules are perfectly 
clear, and then to define it for non-integer s by analytic continuation. This gives 
answers in agreement with the matrix model. 26 ' 27 

Now we show how the multiplicative renormalization A — > /j,q takes place. Al- 
though it seems that eq. (5.3) scales as A s , the free field correlator multiplying it 
is, formally, of order s . For integer s this can be shown through direct calcula- 
tion. Consider, for instance, the tree level correlator of M tachyons with qM < 0, 
9i, ... , qM-i > 0. If we impose the condition s = 0, then qM = 2 — M. The integral 
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for the coefficient of ^(0) can be performed explicitly 27 and gives 

( m - 3)\ l\ r(|«|) ' 

where we have set g\ = 7r 3 /3 2 . If we now send 51, 52, . . . , q m — > 0, while keeping 
Qi = M — 2, we obtain a tachyon correlator with s = m. From eq. (5.5) we 

obtain 

M \ M-l , _ , 

n Tfe) = I^-«+"T(- m )A»0™— 1_ J] ^^^f 1 (5-6) 

\i=m+l I V ; ' i=m+l VIWV 

The physical meaning of the zero of order m is that the "tachyon" decouples at zero 
momentum. This is a sign of the fact that the "tachyon" is actually a massless particle 
in two-dimensional string theory. We will show this more explicitly later on. 

Now we would like to argue that in a cut-off theory s should be replaced by 
(1/1 In jUo|) s - If the cut-off is imposed, then the Liouville zero mode is enclosed in a box 
of size ~ I ln/io|. Thus, the lowest accessible value of the "Liouville energy" 2 + e(q) 
is ~ 1/1 In /iq |, and there is no complete decoupling of T(0). Hence, the amplitudes 
should scale as (A/|ln/zo|) s = (J>q, which is the sought for renormalization of the 
cosmo logical constant. We have identified the physical origin of this renormalization: 
the decoupling of the massless particle at zero momentum. 

The continuation of eq. (5.6) to non-integer s gives the iV-point tachyon ampli- 
tudes with < 0, 91, . . . , qN-l > 0, 27 

This agrees with the matrix model results (4.15)-(4.17) if the matrix model and 
Liouville theory operators are related by momentum dependent normalization factor 

T(q) = f7j^ P (<?) ( 5 - 8 ) 

With this correspondence, we can use the matrix model results to obtain tachyon 
correlation functions to all orders in the genus expansion. This is remarkable because 
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direct Liouville theory calculations of correlation functions have not yet been per- 
formed even for genus 1. For h > 3 we do not know the region of modular integration 
and could not perform a direct Liouville calculation even in principle. From eq. (4.15) 
we find, for instance, 



(T(q)T(k)) = -S(q + k) 



r(i-|g|) n2 
r(M) 



W 91 



i i 2 



,| 24( / 3,) 2 (l ^ 1) ^-^- 1) + - 



(5.9) 

A natural continuation of this formula to q = is to replace lim^o I <? I by 1/| In /i.o | - 
Then we obtain, 

< T (°) T (°)> ~ ^2 E ■ ( 5 - 10 ) 

This is a simple consistency check: differentiation of the path integral with respect 
to A inserts a zero-momentum tachyon. 

To conclude this section, we will describe the only h > calculation that has 
been performed to date. This is the sum over toroidal surfaces with no insertions. 28 
The reason this calculation is relatively straightforward is that s = 0, and the com- 
plicated insertion of T(0) s is absent. Sending s — > in eq. (5.3), and remembering 
the renormalization A — > fiQ, we obtain the free field path integral 

Z x = -iln/xo j d 2 r J [DX][D^)][DbDc]e' So (5.11) 
T 

where r = t\ + iti is the modular parameter of the torus, and T denotes the funda- 
mental region of the modular group: T2 > 0, \t\ > 1, — \ <T\<\. This path integral 
is standard, and we find 

Z^R/V^) = - J In w J d\ {^f) (2n^- 2 \r ] (q)\ 2 )- 1 Z(R/V^,r,f) , (5.12) 

T 

where r\ is the Dedekind function and q = e 2mT . In the above integrand, the first 
term is from the ghost determinant, the second is from the Liouville determinant, and 
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Z(R/y/a', r, r) is the partition function of the scalar field compactified on a circle, 
Z{R Vd,T,T) = ■= ■= , . ., 9 > exp !— L . 5.13 

The double sum is over the soliton winding numbers about the two non-trivial cycles 
of the torus. From eqs. (5.12) and (5.13), we obtain 

Z^R/Vd) = -iln/xo— ^F(R/V^) , 

47tv a' 

^2 / nR 2 \n-mT\ 2 \ ( 5 - 14 ) 

j- z n,m v 7 

The ^-functions cancel out, and we end up only with the contributions of the zero 
modes. This is due to the absence of particles corresponding to the transverse string 
excitations, and reflects the two-dimensional nature of this string theory* Remark- 
ably, the integral in eq. (5.14) can be easily performed using the trick of ref. 30. 
The idea is to trade the m-sum over the winding sectors of the string for a sum over 
many inequivalent fundamental regions which together cover the strip — ^ < t\ < \ 
in the upper half-plane. In 26 dimensions this gives a representation of the string free 
energy in terms of the free energies of all the modes of the string 31 . Similarly, eq. 
(5.14) becomes 

00 

^/^, r~,^ f d 2 r f dr 2 ( nR 2 k 2 \ 
W V50=/ ? + 2y-|pxp(-— j . (5.15) 

The second term is simply the temperature-dependent one-loop free energy for a 
single massless boson in 2 dimensions, expressed in the proper time representation. 
The first term is the one-loop cosmological constant of the massless boson, which has 
been automatically supplied with a 'stringy' ultraviolet cut-off: here the r integral 



★ The isolated transverse states at discrete momenta, to be discussed in section 6, do not affect 
the torus partition function. 29 



27 



is over the fundamental region, not over the strip. Performing the integrals in eq. 
(5.15), we find F(RfVa') = f (1 + a'/R 2 ) and 

Z^R/V^) = "^(^= + ^) In/^o • (5.16) 

As we will show in section 9, this answer agrees with the matrix model result. This 
calculation gives us hope that the Liouville path integral exactly describes sums over 
geometries beyond the tree level. 

6. THE SPECIAL OPERATORS. 

The most striking feature of the tachyon correlation functions is the occurrence 
of the external leg factors ^p^~^ , which contain factorized poles whenever any of 
the momenta approaches an integer value. Although this exact factorization is not 
yet fully understood, it was noted in ref. 20 that the poles are related to the presence 
of other physical states in the theory besides the tachyons (5.2). These states are 

remnants of the transverse excitations of the string and occur only at integer q. 

20,32,33 

The origin of these states can be traced to the pure c = 1 conformal field theory 
where there are special primary fields of the form 34 

Vj, m (dX,d 2 X,...)e 2miX ^ (6.1) 

with conformal weight J 2 . These states form SU(2) multiplets with spin J and 
magnetic number m. The connection with SU(2) becomes apparent once we note 
that the states (6.1) are the full set of primary fields in the compact c = 1 theory 
with the self-dual radius R — 1, where there is a well-known level 1 SU(2) current 
algebra. For each spin J, Vj ; ±j = 1, which gives a tachyon operator. Other mem- 
bers of the multiplet can be constructed by applying raising and lowering operators 
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§ dzexp(±2iX(z)). After coupling to gravity, we can "dress" the special c = 1 
operators to obtain new (1,0) operators of the form 

Vj, m (z) = Vj, m e 2miX e^ J -^(z) . (6.2) 

Recently, it was established that these operators are related to a Wi +OQ algebra. 35 ^ 
Undoubtedly, this algebra plays an important role in determining the properties of 
the theory. 

In the matrix model, the physical (1, 1) operator 

J d 2 zVj, m (z)Vj, m (z) (6.3) 
was, up to normalization, identified with 10 ' 36 

oo 

J dxe 2miX (\ - \ C ) 2J (6.4) 

— oo 

where A c is the coordinate of the quadratic maximum. Using this identification, the 
correlation functions of all operators with m = were calculated in ref. 10. Recently, 
it was shown 37 ' 38 that the operators (6.4) generate a W\ +00 algebra. The connection 
of this algebra with the algebra of special operators in Liouville theory has been 
established by E. Witten. 35 

The simplest new operator is the zero-momentum dilaton 

Vi, Vi,o = dXdX . (6.5) 

Thus, of the full dilaton field of the critical string theory only its zero-momentum 
part remains here. Similarly, the higher J excitations correspond to remnants of the 
higher mass particles. The connection of the special operators to the divergences in 



f A related observation was made by I. Klcbanov and A. Polyakov. 
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the correlation functions is provided by the fusion rules. There are terms ~ l/\z — w\ 
in the fusion rules of the physical operators (6.3). Integration over one operator near 
another can thus produce logarithmic divergences. 

Let us elaborate on this mechanism using, as an example, the tachyon 2-point 
function. 20 The poles in it occur when \q\ is integer, i.e. when tachyon operators 
become members of the SU(2) multiplets of the special operators. The operator 
product expansion of two integer momentum tachyons with n > is 

e inX +{ -2 +H )^ z )e -inX + (-2+\n\)* {Wj + ^_ V|n| _ 1>0 ^ n| _ 1>0 + . . . . 

(6.6) 

Integrating over z near w, we obtain a logarithmic divergence for any integer \q\ > 
due to the appearance of a physical operator in the fusion rule. This explains the 
infinite sequence of poles in the two-point function. 



7. STRING FIELD THEORY FROM THE MATRIX MODEL. 

The success of the Liouville approach shows that the theory of surfaces embed- 
ded in one dimension can also be viewed as critical string theory embedded in two 
dimensions. 7 ' 21 The field 0, which starts out as the conformal factor of the world sheet 
metric, assumes the role of another embedding dimension similar to X. Indeed, the 
path integral (5.1) can be thought of as the sigma model for the bosonic string prop- 
agating in two-dimensional Euclidean target space with the metric G^ u = 5^ u , the 
dilaton background D(<f>) = —40, and the tachyon background T(<f>) = Ae~~ 2 ^. The 
two-dimensionality of the string theory has many important physical consequences. 

First of all, the mass of the particle corresponding to motion of the string in its 
ground state, "the tachyon" is 



Thus, D = c+1 = 2 is the "critical dimension of non-critical string theory" where "the 
tachyon" is exactly massless. Indeed, in the calculation of the Liouville path integral 
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on the torus, we found a massless two-dimensional boson propagating around the 
loop. From the decoupling at zero momentum, observed in the study of the tachyon 
correlation functions, we deduce that the theory is invariant under constant shifts of 
the tachyon field. 

The second expected feature of two- dimensionality is the absence of transverse 
oscillation modes of the string. Indeed, if we could pass to the light-cone gauge, then 
all transverse oscillations would be eliminated and the entire spectrum of the string 
theory would consist of one massless field. In reality, we cannot pick the light-cone 
gauge because of the lack of translation invariance in the ^-coordinate. As a result, 
the spectrum of the string contains the transverse excitations constructed in section 6. 
These states are not full-fledged fields, however, because they only occur at discrete 
momenta \q\ = n. Indeed, the string one-loop calculation of section 5 showed that 
there are no massive transverse fields in the theory. Thus, at the minimal level, we 
should be able to formulate the string field theory with a single massless field, perhaps 
coupled to an infinite number of quantum mechanical degrees of freedom. 

In this section we show that a string field theory with the features anticipated 
above can be derived directly from the matrix quantum mechanics. This further 
strengthens the connection of c = 1 quantum gravity with string theory in D = 2. 
The matrix quantum mechanics directly reduces to the exact string field theory, which 
is the second-quantized field theory of non-relativistic fermions discussed in section 4. 
In this section we will simply recast this formalism so that some of its physical features 
become more transparent. First, we will show that the non-relativistic field theory 
can be expressed in the double-scaling limit as a theory of free quasi-relativistic chiral 
fermions. 10 ' 39 These chiral fermions have the kinetic term that is relativistic to order 
gg t , but receives translationally non-invariant corrections of order g S f. Further, we will 
carry out the conventional bosonization 40 of this fermionic hamiltonian, and obtain 
an interacting theory of massless bosons in D = 2 first derived by Das and Jevicki 
with somewhat different methods. 9 In conclusion, we will perform a few manifestly 
finite calculations in the bosonic formalism to show that it is fully equivalent to the 
original representation in terms of non-relativistic fermions. 
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7.1. Chiral Fermions. 



In order to exhibit our general method, we first consider fermions moving in an 
arbitrary potential U(X). When we later take the double-scaling limit, as expected it 
will depend only on the existence of a quadratic local maximum of U(X). The second 
quantized hamiltonian for a system of free fermions with Planck constant l/f3 is 

h= I dX {w^^ + um ^~ M ^~ N) ] ' (7 ' 2) 

where /ip is the Lagrange multiplier necessary to fix the total number of fermions to 
equal N. As usual, it will be adjusted so as to equal the Fermi level of the N fermion 
system, /ip = e/v- The fermion field has the expansion 

vJ/(A,i) = 5>^(A)e-^ , (7.3) 

i 

where tpi(X) are the single particle wave functions and an are the respective annihila- 
tion operators. The time t is now taken to be Minkowskian. 

All the known exact results of the c = 1 matrix model have been derived from 
this nonrelativistic field theory. An important feature of this theory is that it is two- 
dimensional: the field \&(A, t) depends on t, and also on the eigenvalue coordinate 
A. This is the simplest way to see how the hidden dimension, originating from the 
Liouville mode of quantum gravity, emerges in the matrix model. 

Note that, as (3 — > oo, the single-particle spectrum is approximately linear near 
the Fermi surface, 

^K/iF + ^ + O^j . (7.4) 

This suggests that the behavior of excitations near the Fermi surface is quasi-relativistic 
The spectrum of relativistic fermions in a box is, of course, exactly linear. We will 
therefore find non-relativistic corrections arising from the violations of the linearity of 
e n due to higher orders of the semiclassical expansion. As the Fermi level approaches 
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the maximum, /j, — /j, c — up — > 0, the spacing of the spectrum vanishes, u ~ 7r/| ln/x|. 
Thus, we expect that the size of the box in which the quasi-relativistic fermions are 
enclosed diverges ~ | In /x | . 

Let us introduce new fermionic variables ^ i and 
y(\,t) = ^^=[e- i Pf dX ' v ^ +i ^^ L (X,t) + e^I dX ' v ^^/^ R (X,t)} , (7.5) 

y / 2v(X) 

where v (A) is the velocity of the classical trajectory of a particle at the Fermi level in 
the potential U(\), 

v(X) = ^ = y/2{p. F - U(\)) . (7.6) 

We substitute (7.5) into (7.2) to derive the hamiltonian in terms of and ^r. 
We drop all terms which contain rapidly oscillating exponentials of the form exp [ ± 
2i(3 f X v (x)dx\ , since these give exponentially small terms as (3 ~ — * oo and do 
not contribute to any order of semiclassical perturbation theory. For the same reason 
we can restrict the coordinate A to lie between the turning points of the classical 
motion. For [i > 0, we will consider the case where the fermions are localized on 
one side of the barrier. If we considered the symmetric case of fig. 3, we would have 
found two identical worlds decoupled from each other to all orders of the semiclassical 
expansion. We will also discuss \i < where the Fermi level is above the barrier, and 
the fermions can be found on both sides. 

After some algebra we find 

r/2 

U = 2/3h = [ dr t ^ R dr^ R -t^[d T ^ L + -^(d r ^ L d T ^ L + d r ^ R d T ^ R ) 

{ L PV (7.7) 

5(t/) 2xn 



2w 4 



where v' = dv/dr. Here we see that the natural spatial coordinate, in terms of which 
the fermion has a standard Dirac action to leading order in (3, is r - the classical time 
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of motion at the Fermi level - rather than A. An important feature of the hamiltonian 
is that it does not mix ^ and 



T~t = 7~Ll + Hr . (7.8) 

The only mixing between the different chiralities is through the boundary conditions. 
In order to determine the boundary conditions, consider the leading semiclassical 
expression for \I/(A,0), 

*(A,0) = ^^ n a n +5^^ m 6j n , 

n>0 m<0 

2 f 7T (7-9) 

Ipn = ^== COS (P J d\'v n (\') - -) 

v n = ^/2(e n - U(\)) , 

where we have relabeled on = a\ , i > , ctj = b\, i < 0, and « is the number of the 
fermion energy level starting from the Fermi level \ip. Expanding e n about /ip, we 
find 

*fi=\/l(E e2 ™ r/T «» + £ /T 4) , (7 10) 

n>0 m<0 

^ = \/|( E e- 2mnT/T *n + £ e- 2 ™ /T 4) • 

n>0 m<0 

Thus, and ^ are expressed in terms of a single set of fermionic oscillators, a n 
and b m . Semiclassically, they satisfy the boundary conditions 

*r(t = 0) = * L (r = 0), W fl (r = |) = tt L (r = |) . (7.11) 



These ensure that the fermion number current not flow out of the finite interval, i.e., 



that ^ f (r)7i^(r) = ^t^^ — vanish at the boundary. They also guarantee that 
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tyft and \l/ L are not independent fields and that we are including the correct number of 
degrees of freedom. The issue of dynamics at the boundary is quite subtle, however, 
since the corrections to the relativistic hamiltonian in eq. (7.7) blow up precisely 
at the boundary. The problem is that the semiclassical approximation breaks down 
at the points where v(t) = even as f3 — > 0. One possibility of dealing with this 
breakdown is to carefully regularize the physics at the boundary points. This has 
lead to some success. 41 ' 42 We will suggest another method, which is problem-free 
from the beginning: to approach double-scaling limit with /x < 0. As we will show, 
an equivalent procedure is to work with /i > but to construct the theory in terms 
of ty(p,t), where p is the momentum conjugate to A. 

We have succeeded in mapping the collection of N nonrelativistic fermions, which 
describe the eigenvalues of $, onto an action which, to leading order, is just the two- 
dimensional Dirac action with rather standard bag-like boundary conditions. How- 
ever, the 4 corrections in eq. (7.7) cannot be disregarded in the double scaling limit 
because near the quadratic maximum v(t) = v / 27tsinh(r). In the double scaling limit 
the surviving correction to the relativistic hamiltonian is of order g st = 



dr 



sinh (t) 



|^| 2 + i|^| 2 (l-|coth 2 (r)) 



(7.12) 



where i runs over L and R. This correction does not change the non-interacting 
nature of the fermions, but it does render the fermion propagator non-standard. 

In the double-scaling limit the boundary at T/2 recedes to infinity and is irrele- 
vant. In fact, the correct expression for the hamiltonian would follow had we replaced 
U(X) by — ^A 2 . However, the hamiltonian (7.12) still suffers from the divergence at 
r = 0. In view of this, let us modify the problem to approach the double-scaling 
limit with /i < 0. Since string perturbation theory is in powers of (Pfi)~ 2 , analytic 
continuation to positive /i should not be problematic. Now f 2 (r) = 2|/i| ch 2 r, and 
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we find 



■DO 




dr iW R dr^/ R - iW L d T V L + 



4P\fi\ ch 2 



1 



{d T m\d T m L + d T V R d T * R ) 



— oo 



1 



(l-|th 2 r)(^ L + ^ fl ) 



+ 



8/?|//| ch 2 



(7.13) 



By choosing /i < we have eliminated two problems at once. First, now the hamil- 
tonian has no divergence at r = 0. Second, there is now no relevant boundary at 
all, \&£, and do not mix, and there is no need to impose boundary conditions. 
The fact that the two chiralities do not mix to all orders of the string perturbation 
theory is remarkable. Its impact on the scattering processes will be discussed after 
we introduce the bosonized formalism. 

7.2. The Bosonic Formalism 

In this section we will explicitly confirm the expectation that the space-time 
picture of the c = 1 string theory involves the dynamics of one massless scalar field 
in two dimensions. To this end we will bosonize the fields and following the 
standard bosonization rules for Dirac fermions. 40 A two dimensional free massless 
Dirac fermion is equivalent to a single free massless scalar boson. In our case, although 
the fermions are free, they are not truly relativistic beyond the semiclassical limit. 
This will give rise to interaction terms in the equivalent bosonic field theory. 

To bosonize the system we replace the fermion fields by 40 



where X is a massless two-dimensional periodic scalar field, and P is its canonically 
conjugate momentum. The normal ordering in eq. (7.14) is in terms of the con- 
ventionally defined creation and annihilation operators, which we will utilize for the 



T 




(7.14) 



T 
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explicit calculations of section 6.3. To convert eq. (7.7), we make use of the following 
easily derived expressions 



: *[d T V L - tf R d T * R : 
: V[V L + tf R * R : 
d T *[d T * L + d T tf R d T * R : 



= -:P 2 + (X 
X' 



l\2 



-yft:PX'P + hx'? + ±-X"' 

6 07T 



(7.15) 



Substituting these into eq. (7.7), we find 



T/2 



H :-- 



dr : 



o 
1 



2(3^ 



P l + (X') 



l\2 



71 



1 



1 



X' ( -77 



i\2 



5(v') 
2v 4 



(7.16) 



If we integrate by parts and discard the boundary terms, this reduces to [10] 



H 



T/2 



dr : 



P 2 + (X 



l\2 



n 1 .., •■ I 



(px'p + -(x'Y)- 



2(3^ X \3v 3 



v" (v') 2 



2v 4 



(7.17) 

The boundary conditions obeyed by the field X are determined by those of eq. 
(7.11), which ensure that fermion number not flow out of the r box. As we have shown, 
the current density ^/(t)^i^(t) = ^ r ^r — ^ l&l vanishes at the boundary. Since 
this density is proportional to dtX{t,r), we deduce that X(t,0) and X(t, -|) must 
X satisfies Dirichlet boundary conditions. The constraint on the 



be constant, 



i.e. 



total fermion number requires that X(t,0) - X(t,±) = 0. We are free to choose 
X(t,0) = X(t, -j) = 0. These boundary conditions eliminate all the winding and 
momentum modes of X. As a result, there is no need to worry about the periodic 
nature of X: it acts like an ordinary scalar field in a box with Dirichlet boundary 
conditions. 



37 



The hamiltonian (7.17) and the boundary condition agree with the collective 
field approach. 9 ' 41 Physically, the massless field X describes small fluctuations of the 
Fermi surface. In the double-scaling limit the hamiltonian reduces to (for /i > 0) 



00 





P 2 + (X') 2 ¥ — (px'P + -(X'f) ^= ^X' 

y ! 2/3/1 sh 2 r V r !) UPfiyfF sh 2 r 



(7.18) 

This hamiltonian, as its fermionic counterpart, suffers from a divergence at r = 0. 
This divergence was regularized in ref. 41 with zeta-function techniques. Alterna- 
tively, we may take the double-scaling limit with /i < 0. Then 

oo 3 2 

:H:=- fdr: P 2 + (X') 2 -^-—(PX'P + -(X'f) o ~ 2 ^ -1= X' 

2J I K J 2/5|/i|ch 2 r V 3 V 12(3\fi\ ^ ch 2 r 



(7.19) 

and the sickness at r = has disappeared. The bosonized theory of non-relativistic 
fermions has a remarkable structure. In the double-scaling limit there is a single cubic 
interaction term of order g st , and a tadpole of order g st . Therefore, we can develop 
an expansion of correlation functions in powers of g s t using conventional perturbation 
theory. As we will show, this expansion reproduces the genus expansion of the string 
amplitudes. 



7.3. Scattering Amplitudes. 



Scattering amplitudes of the A-quanta can be related to the Euclidean correlation 
functions in the matrix model. 43 We will exhibit this relation for /j, > 38 (the // < 
case works similarly). The finite boundary operator 0(1, q) from eq. (4.13) can be 
translated into 

oo 

X - J dxe iqx J dre- lX ^ : ^ L + ^ R (r,t) : , (7.20) 
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where A(r) is the classical trajectory at the Fermi level. Upon bosonization, we find 

oo 

0(l,q)~ J dxe iqx J dre- lX ^d T X ~i J dkF{k,l)kX{q,k) (7.21) 



where 



oo 

F(k,l) = J dre- lx ^ cos(kr) , 


X(x,t) = J dxe~ iqx J dksm(kr)X(q,k) . 



(7.22) 



Evaluating F(k, I) with A(r) = ch r, we find 38 

F(k,l) = K ik {ly/2ii) = 2s J( ik7r) ( I -ik(lV^) - hk(lV^)) ■ (7-23) 



For small z, we will replace 



I v (z) - {z/iy^—^ . (7.24) 



In calculating the Euclidean correlation functions, each operator will be connected 
to the rest of the Feynman graph by the propagator l/(q 2 + k 2 ). We will deform 
the /c-integral of eq. (7.21) for each external leg, and pick up the residue of the 
propagator pole. There is a subtlety here: each ^-function has to be split into a 
sum of two /-functions as in eq. (7.23), and the allowed sense of deformation for the 
two is opposite. For I_ ik we pick up the pole at k — i\q\, while for I ik - the pole at 
k — — i\q\. As a result, we obtain the amputated on-shell Euclidean amplitude times 
a factor for each external leg, which for small I is 

''^'''^TTl = -«^)"T(-k|) . (7.25) 
Therefore, the correlation function of puncture operators is 

N \ N 

u=l / j=\ 

where A is the Euclidean continuation of a scattering amplitude of iV X-quanta. The 
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same factor for each external leg appears in the direct fermionic calculations, (4.15)- 
(4.17). The lesson from this rather technical exercise is that the Das-Jevicki field 
theory assembles the correlation functions in a rather remarkable way. Each scattering 
amplitude of the X-quanta has to be multiplied by the external leg factors which arise 
because of the unusual form (7.21) of the external operators. 43 On the other hand, 
these factors containing the poles have an important string theoretic meaning and 
certainly cannot be discarded. Thus, we are still missing the precise interpretation of 
the Das-Jevicki field X in the string theory with Euclidean signature. 

After continuation of string amplitudes to the Minkowski signature, X — > it, the 
connection with the Das-Jevicki field theory is more straightforward. The Euclidean 
external leg factor for the operator T(q) is — ^ q ^ 2 Y{— \q\)/T(\q\). Upon continuation 
\q\ — > — iE, this factor turns into a pure phase 13 



which does not affect any scattering cross-sections. Since the external leg phases are 
unobservable, they can be absorbed in the definition of the vertex operators. Thus, 
we find that the Minkowskian scattering amplitudes of tachyons are given by the 
corresponding amplitudes of the X-quanta. 

In fact, we may consider two kinds of scattering experiments. The usual scattering 
involves colliding right-moving and left-moving wave packets. In this theory such an 
experiment yields trivial results since H = Hl + Hr 43 , and the wave packets pass 
through each other with no influence or time delay. 44 Thus, there is no conventional 
"bulk" scattering, whose rate is finite per unit spatial volume when plane waves are 
being scattered. 

We can, however, consider another kind of scattering. 44 For /x > we may prepare 
a left-moving wave packet incident on the boundary and wait for it to reflect. This 
can be interpreted as the scattering of n left-moving particles into m right-moving 
particles. (Recall that the number of massless particles in one spatial dimension is 




(7.27) 
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a subtle concept that needs to be carefully defined.) For fi < 0, a left- (right- ) 
moving wave packet stays left- (right- ) moving to all orders of perturbation theory in 
g s t. However, now there are two asymptotic regions, r — > ±00, and the wave packet 
undergoes some deformation as it passes through the interaction region near r = 0. 
This deformation can be interpreted as a change of state (and number) of particles. 
It is this "non-bulk" scattering, 44 ' 41 ' 42 whose rate is not proportional to the spatial 
volume, that gives, upon Euclidean continuation, the matrix model correlators. To 
demonstrate this explicitly, we will now calculate tree-level amplitudes for 2 — > 1 and 
2 — * 2 particle scattering. 

We will work with the fi < hamiltonian of eq. (7.19) and show that the calcula- 
tions are manifestly finite. Since the chiralities do not mix, we will consider scattering 
of right-moving particles described by the hamiltonian 



Following ref. 41, we will perform our calculation in the hamiltonian formalism. This 
is advantageous here because bosonization maps the fermionic hamiltonian into the 
normal-ordered bosonic one. Since fermionic calculations are finite, we expect that 
the normal ordering will remove all the divergences in the bosonic approach. In other 
words, the hamiltonian approach provides a definite set of rules on how to handle the 
bosonized theory so that it is perfectly equivalent to the original fermionic theory. 

We introduce the canonical oscillator basis 




00 



1 - f th 2 r 



12/5|/i|v^Fch 2 r 



(P - X>) 



—00 



(7.28) 



OO 



X 



( t)T)= J ^_^ (A; )e^-l^)+at(A;) e ^-l^)), 

—00 



(7.29) 



00 




—00 



such that [a(k), cJ(k')] = 5(k — k'). The hamiltonian assumes the form Hr = H2 + 
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#3 + Hl 



oo 

H 2 = J dkka\k)a{k) , 
o 

00 

H 3 = nA % o\ i / dk\dk 2 dkj,\/ k\k 2 kj, [f(h + k 2 + k s )a(ki)a(k 2 )a(k 3 )- 

24tt/3|/x| J 1 (7.30) 



3f(ki + k 2 - k%) : a(ki)a(k 2 )a) ' (k 3 ) :] + h.c. , 

oo 

Hi = J dkVkg{k){a{k) - a\k)) , 

o 

where 

oo 

ikr 



oo 

1 ^ ixk 



m = J dT dTv e 



— oo 
oo 



9(h) = j 



sh (TTfc/2) ' 

!-| thV ^r_ vr(A; 3 + 2A;) 



(7.31) 



ch 2 r 4 sh (ttA;/2) ' 

oo 

Now we calculate the S'-matrix, 

S = 1 - 27ri5(£i - £/)T . (7.32) 

Each right-moving massless particle has energy equal to momentum, E = k. Consider 
the amplitude for two particles with momenta k\ and k 2 to scatter into a single particle 
of momentum k 3 , 

S(k h k % k 3 ) = -2ni5{E 1 + E 2 - E 3 ) y /E 1 E 2 E 3 < k^H^k^ > (7.33) 

where \k\k 2 >= a) {ki)a) {k 2 )\Q >. Substituting the cubic interaction term, we find 

S(E 1 ,E 2 ] E 3 ) = -g st 5{Ei + E 2 - E 3 )E 1 E 2 E 3 . (7.34) 

The fact that this amplitude describes "non-bulk" scattering is evident from the 
absence of a separate delta-function for momentum conservation. In order to relate 



-e 
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this to the matrix model three-puncture correlation function, we perform Euclidean 
continuation Ej — > i\qj\. Including the external leg factors according to eq. (7.26), 
the result is in agreement with the tree-level contribution to eq. (4.16). 

Now we proceed to the more complicated case of non-forward scattering of par- 
ticles of momenta k\ and k 2 into particles of momenta ks and k 4 . The amplitude is 
given by second-order perturbation theory, 

rrfu u u u\ / tp p p p < hh\H int \i >< i\H int \kik 2 > . , 

T(k h k% fc 3 , fc 4 ) = ^/E 1 E 2 EsE 4 > , (7.35) 

^-^ £4 + t/2 — £>i 

where the sum runs over all the intermediate states i. The s-channel contribution is 
easily seen to be 41 



A CXJ 

2 r 

s^(h, k 2 - h, h) = n Ej j dk 



j- 1 o 



kf 2 (h + k 2 -k) kf 2 (h + k 2 + k) 



k\ + k 2 — k + ie k\ + k 2 + k — ie 

(7.36) 

where the first term arises from the one-particle intermediate state, and the second 
one - from the five-particle intermediate state. Eq. (7.36) can be expressed as an 
integral over k from — oo to oo. The t- and w-channel contributions assume a similar 
form. Summing over the three channels, we obtain 41 

2 4 

S(h,k 2 ;k 3 ,h) = -i^S(E 1 +E 2 -E 3 -E 4 )l[E j (F(p s )+F(p t )+F(p u )) , (7.37) 
where p s = k\ + k 2 , p t = \h - ks\, p u = \ fa - h\, and 



F(p) = / dk 



(p-kf 



sh 2 (7r(p- k)/2) 



k 



p — k + ie sgn(/c) 



(7.38) 



Using 



x — ie x 



V- + m5{x) (7.39) 
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Fig. 4a) The one-loop O {git) Fig. 4b) The tadpole graphs 

contributions to the 1 — > 1 of order g 2 at . 

amplitude. 

we find F(p) = — i^- — Jp. The total amplitude thus becomes 

2 4 

S(E h2 ;E 3A ) = -^5(E 1 + E 2 -E 3 -E 4 )l[E J (E 1 + E 2 + \E 1 -E 3 \ + \E 1 -E 4 \-2t) . 

i=i 

(7.40) 

Upon the Euclidean continuation Ej — > i\qj\, and inclusion of the external leg factors, 
this precisely agrees with the non-relativistic fermion calculation of the 4-puncture 
correlator (4.17). 

Now we give an example of one- loop calculation. Following ref. 41, we calculate 
the 0(g 2 t ) correction to the two-point function (1 — >• 1 amplitude). Like the tree- level 
4-point function, this is given by second-order perturbation theory. The contribution 
to T(k, k) from the one-loop graphs is 

oo 







(247r) 2 J \ k — k\ — k 2 + it k + k\ + k 2 — it 



where the first term is from the 2-particle intermediate state, and the second - from 
the 4-particle intermediate state (fig. 4a). After changing variables to s = k\ + k 2 
and k 2 , and integrating over k 2 , this reduces to 



(24tt) 2 



OO 

f dss z f 2 (k-s) = _gl k2 ( . k3 2k2 8\ 
J k-s + it sgn(s) 48tt V 15 7 



-oo 
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where we evaluated the integral similarly to that of eq. (7.38). This is not the 
complete answer because there are also the diagrams of fig. 4b) arising from the 
tadpole term H\. They contribute equally, giving 

OO OO 

^ 2 f dklkl mm = -M. r dkl jj^L = _A x ^ 2 



(247r)(48vr) J -ki + ie " 384 J sh 2 (nk 1 /2) 48vr 15 



(7.43) 

Adding this to the one-loop contribution (7.42), we find 

T(E, E) = -^-E 1 (iE 3 + IE 1 + l) . (7.44) 

487T 

Upon the Euclidean continuation and inclusion of the external leg factors, this agrees 
with eq. (4.15). 

The agreement with the string one-loop calculation gives us further confidence 
that the bosonic string field theory is finite and exact. This implies that the perturba- 
tion series exhibits the (2h) \ behavior, which is unlike the h\ found in the conventional 
field theory. This appears to be connected with the lack of translational invariance of 
the hamiltonian and of the scattering amplitudes. 41 In general, owing to the expo- 
nential fall-off of the form-factors f(k) and g(k), the hamiltonian (7.30) is completely 
free of ultraviolet divergences. Thus, our innovation of taking the double-scaling limit 
with /i < has rendered the entire perturbative expansion manifestly finite. On the 
other hand, the end results of our calculations do not change under the continuation 
pL — > —fx. This suggests that the divergences for /j, > 0, connected with the fact that 
there f(k) and g(k) did not fall off as k — * oo, 41 are simply a result of using an 
inconvenient formalism. This is also indicated by the fact that all divergences are 
correctly removed by zeta-function regularization. 41 In fact, positive and negative \x 
are related by the transformation that interchanges the classical coordinate A with 
its conjugate momentum p. Thus, the equation 

(4 - = -m ( 7 - 45 ) 
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can be written in terms of p as 

(-3^-3^ = ^- ( 7 - 46 ) 

The lesson is that, while for \x < it is convenient to regard A as the coordinate, for 
H > the natural coordinate is p. If we define the original non-relativistic fermion 
theory in terms of ty(p,t), then the formalism of this section goes through nicely 
for n > 0: there are no boundaries, and the chiralities do not mix. Thus, all the 
calculations above can also be regarded as /i > calculations which are now manifestly 
finite. Using this approach, the tree-level 4-point amplitude was calculated in the 
lagrangian formalism in ref. 42. Another technique, which works efficiently for all the 
tree amplitudes, but is not easily generalized to the loop corrections, was introduced 
in ref. 44. 

Despite the impressive progress in the formulation of the bosonic string field 
theory of the c = 1 quantum gravity, there are many puzzles remaining. In particular, 
how should we interpret the coordinate r and the field X(r)7 Initially, it was argued 
that r is essentially the zero mode of the Liouville field 0. If so, then it is natural 
to think of X as the tachyon field. However, it was shown recently that r is not 
locally related to the scale factor, but is instead a conjugate variable in a complicated 
integral transform. 38 This transform is necessary to turn a string field theory, which is 
highly non-local in 0-space, into a simple local hamiltonian in r-space. Miraculously, 
the matrix model has automatically provided the coordinates in which the theory 
looks the simplest. However, it does not appear possible to interpret X simply as the 
tachyon. Recall that, in addition to the tachyon, the theory has a discrete infinity 
of other non-field degrees of freedom. It appears that the theory has soaked up the 
discrete degrees of freedom together with the tachyon into a single massless field X. 
The form of the matrix model operators for the discrete observables, eq. (6.4), indeed 
shows that they are mixed in the X-field together with the tachyon. Clearly, we need 
to attain a much better understanding of the precise string theoretic meaning of the 
bosonized field theory of the matrix model. 
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8. COMPACT TARGET SPACE AND DUALITY. 



In this section we consider the discretized formulation of the sum over surfaces 
embedded in a circle of radius i?. 6 We will adopt the basic definition (2.5), with the 
new condition 

G(x) = G(x + 2ttR) (8.1) 

that ensures the periodicity around the circle. In this model we encounter a new set 
of physical issues related to the independence of the sum. The c = 1 conformal field 
theory is symmetric under the transformation 45 

^ a' voU . _ 

R ^r> Qst^Qst^- ■ (8.2) 

The change of g s t, equivalent to a constant shift of the dilaton background, is necessary 
to preserve the invariance to all orders of the genus expansion. A symmetry of the 
matter system, such as (8.2), is expected to survive the coupling to gravity. In 
this section we will find, however, that in the discretized formulation the duality is 
generally broken, and can only be reinstated by a careful fine-tuning. This breaking 
of duality is due not to the coupling to gravity, but to the introduction of the explicit 
lattice cut-off. Indeed, a generic lattice formulation of the compact c = 1 model on a 
fixed geometry has no dual symmetry. The phase transition that separates the small 
R (high temperature) phase from the large R (low temperature) phase is well-known 
in statistical mechanics as the Kosterlitz-Thouless transition. 46 Physically, it is due 
to condensation of vortices, configurations that are ignored in the "naive" continuum 
limit, but are included in the continuum limit of a generic lattice theory. The vortices 
are irrelevant for R > Rkt, but they condense and change the behavior of the theory 
for R < Rkt- 

In this section we will demonstrate the lack of duality of the discretized sum over 
surfaces. We will also show that, if the lattice sum is modified to exclude the vortices, 
then the partition function is explicitly dual. These general results will be of interest 
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when, in the next section, we consider the specific discretized sum generated by the 
matrix model. We will be able to isolate the effects of vortices in the matrix model, 
and will calculate the dual partition function in the vortex-free continuum limit. 

In the conformal field theory the target space duality (8.2) is generated by the 
dual transformation on the world sheet 47 

d a X^e ab d b X . (8.3) 

Its lattice analogue is the transformation from the lattice A to the dual lattice A. 
Recall that, when we sum over lattices A, we define the target space variables X{ at the 
centers of the faces of A, i.e. at the vertices of A. We will carry out a transformation, 
after which the new variables, pj, are defined at the vertices of A. 

First, to each directed link < ij > of A associate the link < I J > of A which 
intersects it, directed so that the cross product ij x I J points out of the surface (fig. 
1). Now we define Dij = Djj = x% — xj. Let us change variables in the integral (2.5) 
so that the integral is over Djj instead of X{. There are E links but, apart from the 
zero mode, there are only V — 1 independent x's. Therefore, the Djfs are not all 
independent but must satisfy F — 1 + 2h constraints. The constraint associated with 
each face is that Yltij) Dij = 0, where the sum runs over the directed boundary of 
a face of A. This is equivalent to the condition that the sum of the D/j's emerging 
from each vertex of A must vanish, Ylj^u = 0. Similarly, for each independent 
non-contractable loop on A, of which there are 2h, we find Ylioop-^v = 0- ^ n terms 
of A this means that tfj^u = 0, where the symbol ejj is non-zero only if the 

link IJ intersects the specially chosen non-contractable loop a on A. We direct the 

— * 

loop a, and define efj — 1 if I J x a points into the surface, and —1 if it points out. 

The theorem of Euler V — E + F = 2 — 2h, insures that the net number of 
independent variables remains unchanged. Introducing Lagrange multipliers pi for 

* Here V, E and F are, respectively, the numbers of vertices, edges and faces of A, and h is the 
genus. 
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each face of A and l a for each non-contractable loop, the integral over the variables 
Xi, for each discretization of the surface, can be replaced with 



J J\dD IJ G{D I j) J Y[d Pl exp 



{U) 

2h 



U 



2h 



Y[ dl a exp 



a=l 



il a J2t a ijDn 

(IJ) 



I n n di - n - ^ + 

J I a=l (IJ) 

(8.4) 

where G is the Fourier transform of G. We see that, after the transformation (8.4), 
the Lagrange multipliers pi assume the role of new integration variables residing at 
the vertices of A. On surfaces of genus > there are 2h additional /-integrations. 
Geometrically speaking, one introduces a cut on the surface along each canonically 
chosen non-contractable loop, so that the values of p undergo a discontinuity of l a 
across the a-th cut, and subsequently integrates over l a . On a sphere these additional 
integrations do not arise, and one simply includes a factor G{pj — pj) for each link 
of the lattice A. 

The discussion above literally applies to the string on a real line. When G is 
defined on a circle of radius R, then the dual variables l a and pj assume discrete 
values The integrals over pj and l a in eq. (8.4) are replaced by discrete sums, 
and the partition function, when expressed in terms of the dual variables, bears little 
resemblance to the original expression. In fact, the new variables pi describe the 
embedding of the world sheet in the discretized real line with a lattice spacing 1/R. 
On a sphere, where l a do not arise, the resulting lattice sum is precisely what is 
needed to describe string theory with such an embedding. For higher genus, the few 
extra variables l a spoil the precise equivalence, but this does not alter the essential 
"bulk" physics. 

We have found that the transformation to the dual lattice, unlike the transforma- 
tion (8.3) in the "naive" continuum limit, does not prove R — > a' / R duality. Instead, 
it establishes that string theory embedded in a compact dimension of radius R is 
essentially equivalent to string theory embedded in a discretized dimension of lattice 
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spacing 1/R. To conclude that R is equivalent to 1/R is just as counterintuitive as 
to argue that theories with very small lattice spacing in target space have the same 
physics as theories with enormous lattice spacing. In section 11 we will directly study 
the theory with the discrete target space. We will show that, if the lattice spacing 
is less than a critical value, then the target space lattice is "smoothed out", but if 
it exceeds the critical value, then the theory drastically changes its properties. The 
two phases are separated by a Kosterlitz-Thouless phase transition. Equivalently, in 
the circle embedding the massless c = 1 phase for R > Rkt is separated by the K-T 
transition from the disordered (massive) c = phase for R < Rkt- The mechanism 
of this phase transition involves deconfinement of vortices. 46 

Consider a fixed geometry g^ V) and assume that in the continuum limit the action 
reduces to eq. (2.2). A vortex of winding number n located at the origin is described 
by the configuration 

X{9) = nR6 (8.5) 

where 6 is the azimuthal angle. This configuration is singular at the origin: the values 
of X have a branch cut in the continuum limit. On a lattice, however, there is no 
singularity, and the vortex configurations are typically included in the statistical sum. 
If we introduce lattice spacing ^/Ji, then the action of a vortex is S n = n 2 R 2 \ In /x|/4a', 
and it seems that the vortices are suppressed in the continuum limit. We will now 
show that this expectation is only true for large enough R. Let us consider the 
dynamics of elementary vortices with n — ±1. Although each one is suppressed by 
the action, it has a large entropy: there are ~ l//i places on the surface where it 
can be found. Thus, the contribution of each vortex or antivortex to the partition 
function is of the order 

Ie" Sl = ^ 2 /4a')-l . (8.6) 

It follows that, for R > 2v / c7, the vortices are irrelevant in the continuum limit. On 
the other hand, for R < 2\/c/ they dominate the partition function, invalidating the 
"naive" continuum limit assumed in the conformal field theory. 46 In fact, in this 
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phase the proliferation of vortices disorders X and makes its correlations short range. 
The critical properties of this phase are those of pure gravity. Is there any way to 
salvage the "naive" continuum limit? The answer is yes, but only at the expense of 
fine tuning the model, so that the vortices do not appear even on the lattice. Below 
we find an explicit example of a model without vortices, and show that it does exhibit 
the R-*l/R duality. 

First, we have to give a clear definition of a vortex on a lattice. Loosely speaking, 
a face / of lattice A contains v units of vortex number if, as we follow the boundary of 
/, the coordinate x wraps around the target space circle v times. This does not quite 
define the vortex number because x is only known at a few discrete points along the 
boundary and cannot be followed continuously. In order to define vortex number on 
a lattice, it is convenient to adopt the Villain link factor 48 

oo 

G(x) = e-^ x+2nmR f (8.7) 

m=— oo 

where the sum over m renders G periodic under x x+2ttR. Now, eq. (2.5) becomes 

V 27vR oo 

^o,-) = E^ 1) E- Area n /^=II E e-^-^ 2 —*) 2 . (8.8) 

h A i=l l (ij) 

In this model the number of vortices inside a face / can be defined as 

V j = ^rriij , (8.9) 
dl 

where dl is the directed boundary of /. After the dual transformation, one obtains a 
discrete sum with 

G{pi -pj + l a e a u ) = e~^-^+ 1 ^ 2 

m ? n a ■ (8-10) 

PI= R> L= R 

This sum is closely connected to the sum over all possible vortex numbers inside each 
face of A. 48 
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We shall now change the definition of the partition function (8.8) to exclude the 
vortex configurations. It is helpful to think of link gauge fields. The vortex 

numbers vj defined in eq. (8.9) are then analogous to field strengths. If there are no 
vortices, the field strength is zero everywhere. However, we still need to sum over all 
possible windings around non-contractable loops on the lattice A, I a = Ylioop a m ij- 
Thus, the space of rriij we need to sum over is 

rriij = efjlA + mi - rrij, (8.11) 

where mi range over all integers and play the role of gauge transformations, ef- is 
defined above and is non-zero only for the links (ij) intersecting the specially chosen 
non-trivial loop A on the lattice A. Thus, I a is the winding number for the non-trivial 
cycle a on A which intersects the loop A. Summing over rriij from eq. (8.11) only, 
the partition function becomes 

V-l °° oo 

Z( 9o ,.) = 2,Rj29 2 o {h - 1) j:- V ll J J| £ Y[ e -^-^™^\ 

h A i=l_ 00 Ia=— oo (ij) 

(8.12) 

where the factor of 2irR arises from integration over the zero mode of x. After 
transforming to the dual lattice we obtain 




(pi-pj+Lelj/R) 2 



(8.13) 

Thus, after elimination of the vortices, the transformation to the dual lattice clearly 
exhibits duality under 




(8.14) 



Remarkably, the duality is manifest even before the continuum limit is taken. In 
the next section we consider the matrix model for random surfaces embedded in a 
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circle. We will find a clear separation between the vortex contributions, and the 
"naive" vortex-free continuum sum. This will allow us to calculate explicitly the dual 
partition function of the vortex-free model. 



9. MATRIX QUANTUM MECHANICS 
AND THE CIRCLE EMBEDDING. 

The Euclidean matrix quantum mechanics is easily modified to simulate the sum 
over discretized random surfaces embedded in a circle of radius R: we simply define 
the matrix variable 3>(x) on a circle of radius R, i.e. $(rr + 2nR) = The path 

integral becomes 6 



2tvR , g > 

Z= J D N2 $(x)exv -p J dsTrm^j + £/($) 



(9.1) 



The periodic one-dimensional propagator 

00 

G(xi-Xj)= e-\ Xi - x ^ mR \ (9.2) 

m=— oo 

gives the weight for each link in the random surface interpretation. In terms of the 
hamiltonian (3.10), eq. (9.1) is simply a path integral representation for the partition 
function 

Z R = Tre- 27rR P H , (9.3) 

so that 2nR plays the role of the inverse temperature. The problem of calculating 
the finite temperature partition function seems to be drastically more complicated 
than the zero temperature problem, where only the ground state energy was relevant. 
Now we have to know all the energies and degeneracies of states in arbitrarily high 
representations of SU(N). We also have to explain the sudden enormous jump in 
the number of degrees of freedom as we slightly increase the temperature from zero. 
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In view of the discussion in the last section, the reader should not be surprised if we 
claim that the new degrees of freedom, incorporated in the non-trivial representations 
of SU(N), are due to the K-T vortices. Since the vortices are dynamically suppressed 
for R > Rkt-, the problem of calculating Zr is not as complicated as it may seem. 

Let us begin by calculating the contribution to Zr from the wave functions in 
the trivial representation of SU(N). Since the singlet spectrum is that of N non- 
interacting fermions moving in the potential U, its contribution to the partition func- 
tion can be calculated with the standard methods of statistical mechanics. 

Instead of working with a fixed number of fermions N, we will take the well-known 
route of introducing a chemical potential pp adjusted so that 



_2 _ 



oo 

N 



— ■ / p(e) „ [, , rde . (9.4) 

J rw 14- e 2TvRf3(e-fj, F ) v > 



In the thermodynamic limit, N — > oo, the free energy satisfies 

dF 



SN ^ ■ (9 ' 5 > 



Shifting the variables, \i = ji c — /ip, e = fx c ~ e ; we can write the equations for the 
singular part of F as 



7rp(yLi) = 7r J dep{e 



dp I J r w dp 1 4- e 2irR P0*-e) 

— oo (.9.DJ 

OF 1 
-k~t = -PP , 

OA 71 



to emphasize resemblance with eqs. (3.18) and (3.19) in the R = oo case. Indeed, 
we will calculate F(A) using the same method. The change introduced by a finite R 
only affects the first of the equations, determining p(A). We can think of p(p) as the 
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temperature- modified density of states. Differentiating it, we find 

00 

l - 9 -l = [ de l - 9 -l ^ (9 7) 

Pdn J pde 2cosh 2 [7iRP(fx-e)] ' K ' J 

— oo 

Now, substituting the integral representation (3.26) and performing the integral over 
e, we find 



oo 



1 dp _ 1 T ^ f T/2 T/2R 



-Im / dTe" ipf11 - 1 — ; (Q 8) 

(3dfi Ti J sinh(T/2) sinh(T/2i2) ' 1 J 

o 

up to terms 0(e^ /3 ^) and 0(e^ l3 ^ 2nR ) which are invisible in the large-/? asymptotic 
expansion. Thus, this integral representation should only be regarded as the generator 
of the correct expansion in powers of g s t- Integrating this equation and fixing the 
integration constant to agree with the WKB expansion, we find 

oo 

— = Re/-e-** t/2(3fi t/2(3fiR (99) 
<9/i J t sinht/2/3/i sinh(t/2/3/x J R) ' 1 ' ' 

This relation has a remarkable duality symmetry under 

R^ 4' p^Rp . (9.10) 

R 

This is precisely the kind of duality expected in the vortex-free continuum limit. 
This strongly suggests that, in discarding the contributions of the non-singlet states 
to Zr, we have suppressed the vortices. Later on we will offer additional arguments 
to support this claim. 

To demonstrate the duality in the genus expansion, consider the asymptotic ex- 
pansion of eq. (9.9) 

-2m 



OA 

d/j, 



ln/i+^(2/W^) frn(R) 
m=l 

f (z>\ — (0m m\^l9 2fc 91 |o2(m-fe) PI [jj^j \ B 2{m-k)\ pm _ 2 fc 

f m (R)- (2m -1)\^\2 -2||2< - 2| _ _i2 



(9.11) 



fc=0 

Note that the functions /^(-R) are manifestly dual: /i = l(-R + l/i2), /2 = ^7)(7-R 2 + 
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10 + 7R 2 ), and so forth. Solving for //(A) and integrating eq. (9.6) we find the sum 
over connected surfaces Z = —2nRf3F: 



(2(3fi VR) 2 In u - 2h(R) lnu + V f ^ +l{R ], (2(3^y/Ry 

m (2m + 1 

m=l 



2m 



(9.12) 

Comparing this with eq. (3.30), valid in the case of infinite radius, we find that the 
coupling constant l/((3/io) has been replaced by 

»°» {R) = j^m • (9 - 13) 

and that the contribution of each genus h contains the function fh{R). Since both 
g e ff(R) and fh(R) are invariant under the dual transformation (9.10), we have con- 
firmed that the genus expansion of the sum over surfaces is manifestly dual. As in the 
critical string theory, the effective coupling constant depends on the radius. Indeed, 
the transformation R — > 1/R keeping /3«o fixed is not a symmetry, as it interchanges 
a weakly and a strongly coupled theory. 

From eq. (9.12) we find that the sum over toroidal surfaces is —-^(R+l/R) In «o- 
This agrees with the direct calculation in the continuum, eq. (5.16). The extra 
factor of 2 in the matrix model is due to the doubling of free energy for models with 
symmetric potentials. To get the basic sum over triangulated surfaces, we have to 
divide eq. (9.12) by 2, which then gives perfect agreement with Liouville theory. 
This provides additional evidence that the singlet free energy in the matrix model 
evaluates the sum over surfaces in the vortex- free continuum limit. To see this more 
explicitly, we have to consider the non-singlet corrections to the free energy and show 
that, like the vortices, they are irrelevant for R > Rkt- 

This problem was discussed in ref. 12, where it was argued that the total partition 
function can be factorized as 

where 8E n is the energy gap between the ground state and the lowest state in the nth 
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representation, and D n is the degeneracy factor. Let us discuss the leading correction 
to the free energy coming from the adjoint representation. Although the degeneracy 
factor D a( jj = N 2 — 1 diverges as N — > oo, an estimate of ref. 12 shows that the 
energy gap also diverges in the continuum limit, 

(35E adj = c\\n^\ . (9.15) 

As a result, the correction 

D adj e~ 2nmE ^ -~ N 2 fi 27vRc (9.16) 

is negligible for R > ^ in the continuum limit where N/j, is kept constant. Taking 
the logarithm of eq. (9.14), we find that the total free energy is 

F = F singlet + O (n 2 ^ Rc ) + ... (9.17) 

Thus, on the one hand, the higher representations are enhanced by enormous degen- 
eracy factors; on the other hand, they are suppressed by energy gaps which diverge 
logarithmically in the continuum limit. This struggle of entropy with energy is pre- 
cisely of the same type as occurs in the physics of the Kosterlitz-Thouless vortices, as 
discussed above. In fact, we will now argue that the leading non-singlet correction in 
eq. (9.17) is associated with a single vortex-antivortex pair. We perturb the action 
of the continuum theory with the operator which creates elementary vortices and 
anti- vortices, 

O v = J d 2 a^cos[^(X L -X R )} (9.18) 

where Xl r are the chiral components of the scalar field X. In string theory language, 
this is the sum of vertex operators which create states of winding number 1 and —1. 
The insertion of such an operator on a surface creates an endpoint of a cut in the 

- t-,2 

values of X. The conformal weight of O v is h — h — j^i] therefore, O v becomes 
relevant for R < 2\fa' ', causing a phase transition at Rrt = 2\/o/ 49 . The phase 
transition is connected with the instability towards creation of cuts in the values of 
X, i.e., with vortex condensation. 
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This picture of the Kosterlitz-Thouless phase transition can be easily adapted 
for coupling to two-dimensional quantum gravity. The sum over genus zero surfaces 
coupled to the periodic scalar reduces to the path integral 

Z = N 2 J[DX][D<j>]exp(-S L - S - k J d 2 a y/§e^ cos[^(X L - Xr)]^ (9.19) 

where Sl is the Liouville gravity action for c = 1, and (3 = — 2 + -j= so as to make the 
conformal weight of the perturbing operator (1,1). The gravitational dimension of 
the perturbation is then — 1. For small values of the cut-off and for R > 2\fa' \ 

^ = -f, 2 R Inp + O (k V 7 ^ 7 ) + O (k A ^ 2R l^)- 2 ) + . . . , (9.20) 

which demonstrates that the Kosterlitz-Thouless transition occurs at R c = 2y/a/, the 
same value as in a fixed gravitational background. The term 0(k 2 ) in Z originates 
from the 2-point function of dressed winding mode vertex operators, i. e. from surfaces 
with one cut. This leading correction has the same form as eq. (9.17) found in 
the context of matrix quantum mechanics. The further corrections, coming from 
additional vortex-antivortex pairs, are even more suppressed for R > Rkt m the 
continuum limit. 

For R < Rkt, however, the expansion in k diverges badly, indicating an instabil- 
ity with respect to condensation of vortices. In this phase we expect the field X to 
become massive, and the continuum limit of the theory should behave as pure gravity. 
In the matrix model this is easy to show for very small R. The argument is standard: 
at very high temperature a d + 1-dimensional theory reduces to <i-dimensional theory 
because the compact Euclidean dimension of length 2nR = 1/T becomes so small that 
variation of fields along it can be ignored. Thus, for small R, the matrix quantum 
mechanics (9.1) reduces to the integral over a matrix (zero-dimensional theory) 50 



/ 



D N $exp[—27ri?/3Tr [/(<£»)] , (9.21) 
well-known to describe pure gravity. 
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We have accumulated a considerable amount of evidence that the non-singlet 
wave functions of the matrix quantum mechanics implement the physical effects of 
the Kosterlitz-Thouless vortices. For R > Rkt their corrections to the partition 
function are negligible, and the singlet free energy (9.12) gives the sum over surfaces. 
For R < Rkt the non-singlets dominate the partition function, and the matrix model 
no longer describes c = 1 conformal field theory coupled to gravity. In this phase the 
singlet free energy gives the sum over surfaces in the vortex-free "naive" continuum 
limit. As we showed in section 8, this continuum limit can be achieved at the expense 
of a careful fine tuning of the lattice model. In the matrix model the equivalent 
operation turns out to be simple: just discard the non-singlet states. We conclude, 
therefore, that for any R the singlet sector of the matrix quantum mechanics describes 
c = 1 conformal field theory coupled to gravity. This remarkable fact means that, like 
the R = oo theory, the compact theory is exactly soluble in terms of free fermions. In 
the next section, we will use this to calculate compact correlation functions exactly 
to all orders in the genus expansion. 

10. CORRELATION FUNCTIONS FOR FINITE RADIUS. 

In this section we evaluate the correlation functions of the compact c = 1 model 
coupled to gravity using the formalism of free non-relativistic fermions at finite 
temperature. 13 We have already used this approach to find the sum over surfaces 
with no insertions, eq. (9.12). The calculation of correlation functions has a very 
simple relation to the zero-temperature (R = oo) calculation described in section 4. 
We again consider general operators of type (4.7), which reduce to fermion bilinears. 
Thus, the goal is to calculate the generating function G(qi, X\; . . . ; q n , X n ) of eq. (4.9) 
in the compact case. 

The free fermions are now described by the thermal Euclidean second- quantized 
action 

oo 2-kR 

s =I dx J d ^(-i + ^ + ^ + M<P, (io.i) 

-oo 
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and the fermion field satisfies the antiperiodic boundary conditions tp(X,x + 2nR) = 
—^(X,x). Compactifying the Euclidean time is a standard device for describing the 
theory at a finite temperature T = l/(2nR). The thermal vacuum satisfies 



((3li\a\{v)a e ,{v)W) = <WM = 1^— (10.2) 



where f(y) is the Fermi function and \x is the Fermi level. The thermal Euclidean 
Green function may be obtained by replacing 6{y — fl/i) by f(v) in the T = Green 
function of eq. (4.10), 



S E (x 1 ,X 1 ;x 2 ,X 2 ) = J dve-^-^ Ax {6(Ax)f(v)-6(-Ax)(l-f(v))}x 

( ^ 10.3 

sgn(u> n )oo y ' 



o 



The only change compared to the R = oo formula (4.10) is that the allowed fre- 
quencies become discrete, u n = (n + \)/R-, where n are integers. This ensures the 
antiperiodicity of Se around the compact direction. 

As in the zero-temperature case, the calculation of eq. (4.9) reduces to a sum 
over one-loop diagrams. Since the frequencies are now discrete, the single integration 
over the loop momentum is replaced by a sum. This is the only change introduced by 
the finite temperature. As we show below, its effects are easy to take into account. 

First of all, each external momentum is quantized, qi — ni/R. As a result, 
the momentum conserving delta function 5(^2 qi) occurring in the non-compact case 
is replaced by R times the Kroenecker function S^ qi - Consider a diagram which 
corresponds to a given ordering of the external momentum insertions, which we denote 



GO 



as 91, 92, • • • , In- Factoring out the delta-function, the contribution of this diagram is 

sgn(co n )oo sgn(ui n +qi)oo 



x 



(Ai|e 2 ' tQ2iJ |A 2 )... y do; n e- Q "( a '" + S"" 1 ft)+^«n^ n _ 1 | e 2ia„ J ff| Ari ^ 



(10.4) 

It is necessary to break the sum over oj n into parts where the sgn functions in the 
exponents are constant. Thus, we form Pk = — Yli Ik an d put the p's in the increasing 
order. If p and p ! are two consecutive p's, then the summation that needs to be 
performed is, defining £ = a i, 



P'- 1 



r ^ , e " e sh • (ia5) 

This factor multiplies the integrand of an n-fold a^-integral. Comparing this with 
the R = oo case, we find that there the loop momentum integration between p and 
p' gives the factor (e~ p ^ — e~ p ^)/£ multiplying the same n-fold integrand. It follows 
that, in every single contribution to the final answer, the compactification of the 
target space simply inserts the extra factor ^{^j2R) * n ^° ^ e integrand. Thus, from 
eq. (4.12) we arrive at 



ds 



l 



oo eioo 

°" eE ™ -oo 
e„_ioo 

/ cfan-xe-^--- "^(A^le^A^) . . . (Kin)^-^ « )H \\,(i)) 



(10.6) 

This formula allows us to calculate correlation functions of any set of operators of type 
(4.7) in the compact c = 1 model. The integral representation of any such correlator 
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is obtained from that in the non-compact case by insertion of the factor ^^j^R) 
the integrand, and by replacement S(Y2Qi) ^J2<n' ^ s P ec ^ a ^ case °f this rme ^ s 
evident in passing from the R = oo equation (3.16) to the finite R equation (9.8). 

In general, comparing eqs. (10.6) and (4.12), we find that, if 

qi, . . .,q n ;n,(3) 



then 



(Ci(gi) . . . O n (q n )) R = B£<£ q .F R (qi, . . . , q n ; fi, (3) , 

l d 

FR(qi, ...,q n ;fi,P) = —j^^F(q h ...,q n ;V, P) 



(10.7) 



This remarkable relation exists because the entire dependence on /x is through the 
factor e l @^ in the integrands of eqs. (4.12) and (10.6). Therefore, a function /(£) 
can be introduced into the integrand simply by acting on the integral with operator 

Since /i has the interpretation of the renormalized cosmological constant, in each 
sum over surfaces of fixed area the compactification simply introduces the non- 
perturbative factor g^ppr^g) ■ The remarkable feature of this modification, and, 
equivalently, of the operator in eq. (10.7), is that they are entirely independent of 
which operators are inserted into the surface. To demonstrate the power of eq. (10.7), 
we will use the non-compact correlators of eqs. (4. 15)- (4. 17) to develop the corre- 
sponding genus expansions in the compact case with virtually no extra work. To this 
end, we expand the operator in eq. (10.7) as 



Substituting this expansion into eq. (10.7), we determine the finite R genus g corre- 
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lators in terms of the R = oo genus < g correlators, 
F 9 =\ qi ^) = F^\ qi -ri , 

Ft 1 fan) = F 9=1 to^) + ^R-2^ F9=Q ton) > 

Fr\^) = + ^^\q„) + , etc . 

(10.9) 

In particular, apart from the discretization of momenta and the change in the overall 
(^-function, all spherical correlation functions are not affected by the compactification. 
This was, of course, expected in the Liouville approach. The advantage of the matrix 
model formalism is that the entire genus expansion of any correlator can be calculated. 
Surprisingly, as we have demonstrated, the compact case is hardly any more difficult 
than the non-compact case. From the point of view of the Liouville approach, this 
conclusion is quite remarkable. There the finite R calculations have the additional 
complication of having to perform 2h sums over the winding numbers around the non- 
trivial cycles of a genus h surface. Thus, eqs. (10.7) and (10.9), which are completely 
general, pose a new challenge for the Liouville approach. 

Below we give genus expansions for some correlation functions of puncture oper- 
ators P(q). For the two-point function up to genus three we get 

i (VR^r 2 



(P(q)P(-q)) = -R\T(l-\q\)YiJL 



q\ 24 

4 3 



-(kl-i)x 



W-M-D-^^IItol-Ox 

R 2 (3q 4 - 10\qf - 5q 2 + 12\q\ + 7) - 10(g 2 - \q\ - 1) + 

" o{#W - 63|g| 5 + 42g 4 + 217|g| 3 - 205|g| 



93) 



2903040 

r=I 

1 47 9^? 1 

- 21i?(3g 4 - 10|g| 3 - 5q 2 + 12\q\ + 7) + — (q 2 - \ q \ - 1) - — K + . . . 

"(10.10) 

We see that the external leg factors are the same as in the non-compact case. The 
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relation between the tachyon operators T(q) and P(q) is also the same, eq. (5.8). For 
the same reason as in the non-compact case, eq. (10.10) is, strictly speaking, valid 
only for non-integer q, where we may replace fj, by fj,Q. As in the sum over surfaces, we 
find that the effective string coupling constant is g e ff ~ • The -R-dependence 

of the correlator at each genus is not dual: the duality is broken by insertions of 
momentum into the world sheet. In the limit q\ — > the duality is restored, and the 
tachyon correlators reduce to derivatives of the sum over surfaces Z(ixq,R), 



(T(O)T(O)) 



d 2 
dA 2 



(10.11) 



The three-point function for the case ?i, ?2 > is, up to genus two, 



(T( q i)T(q 2 )T(q 3 )) = R5 Eqi ]J[% T ^^ 1 



r(|*|) 



^ )_2 (k3| - l)(\Qs\ - 2) \R(Ql + £ ~ ks\ ~ 1) - ±\ + Ildftl - Ox 

r=l 



21 " " vri '" ~ | " ' ! _/ i? 

^ 2 ( %1 + <£) + 10 ?i 2 ^ 2 - 10(g 2 + ql)\qi\ - 5(gi - q 2 f + 12\q 3 \ + 7 



-10(q 2 1 +q 2 2 -\q 3 \-l) + 



7 \ (VRf3fi)- A 



+ ... 



R 2 J 5760 

(10.12) 

In the case of the four-point function we must consider two independent kinematic 
configurations which will give amplitudes that are not related by analytic continuation. 19 The 
simplest case is when qi,q2,qs > 0, 
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(T( qi )T(q 2 )T(q 3 )T(q 4 )) = -R5^ qz J] 



i=l 



r(i-|gt 
r(|«l) 



H 3 j (/fy) 



-2 



(l?4|-l) 



24 



n(|9 4 | - r) ji?(g 2 + q\ + g| - |g 4 | - 1) - ^} 

r=l ^ ' 



/ — ^ 

- Ild^l - r)^ 2 ( 3(^ + g 2 4 + 4) + 10(g 2 g 2 + q\q\ + g 2 g 2 ) 
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,2.2, 2\ r/„2 



- I0|g 4 |(gi + ?2 + - %if + d + d) + 10 (^ig2 + gigs + g2gs) + I2|g 4 | + 7 



7 



- I0(gf + g| + gf - |g 4 | - 1) + j + • • • 

(10.13) 

The other kinematic configuration is gi,g2 > 0, g3,g 4 < with min(\qi\) = q\, 
max(\qi\) = q 2 . This leads to a much more complicated formula, so we just give the 
answer out to genus one 

4 

(T( qi )T(q 2 )T(q 3 )T(q 4 )) = -R8^ qi \[ 



i=l 



r(l - |gj|) m\ g 2 - 1 

jJL 2 ' 



r(|ft|) 



i - 



24 



x ^R{S a (S - 2q 2 - 6) + g 2 (S 2 + g 2 + q\) + 13^g 2 - 8S<£ - 5g 2 (^ + qf) 
+ 6S(S - 2q 2 ) + 6(g| + q\) + 10g 2 2 + S - 2q 2 - 6^ - ^(5 - 2)(S - 3)| + . . . 



2 , „2^ 



(10.14) 



where we have defined S — h^2\qi 



Using the matrix model, we have directly calculated the correlation functions of 
the tachyon operators for the momentum modes, 



T{q) = i / Sa^^\ 



;io.i5) 



where X and X are the holomorphic and anti-holomorphic parts of the X-field. In 
string theory we also have the winding operators 



10.16) 



G5 



where q is quantized as nR. Although these operators are hard to introduce in the 
matrix model directly, their correlators are obtained from those of the momentum 
operators through the dual transformation 



X + X -> X - X , 



13 -> PR. 



(10.17) 



Carrying out the dual transformation (10.17) on eqns. (10.10)-(10.14), we find 



f(q)f(-q))=R 



r(kl) . 

{VRP^r 2 

24 



/i 



kil 



1 

\q\ 



1) <! y,(<f 



1 



- \Q\ 



1)-R 



(10.18) 



etc. In general, if 



(T( gi ) . . . T(q n )) R=00 = 5(J2 <li)F{qi, . . . , q n ] H, P) , then 



f(qi) . ..f(q n 



R 



sin 



2/3 9^ 



F(q h ...,q n ;fj,,/3R) 



(10.19) 



This formula was obtained by carrying out the dual transformation on the similar 
formula (10.7) for the momentum states. Since we do not have a direct matrix model 
calculation of the winding operator correlations, we do not yet know the expressions 
for correlators where both the momentum and the winding operators are present. 

Is there any real time interpretation of the exact correlation functions we have 
calculated? In the R = oo case we found that, after the continuation to Minkowski 
signature, the external leg factors turn into pure phases, and the correlation functions 
become scattering amplitudes of the X-quanta of the Das-Jevicki field theory. We can 
perform the same continuation, \qj\ — > —iEj, for the finite R correlation functions, 
expecting to obtain the same field theory at a finite temperature T = l/(2nR). The 
external leg factors again turn into phases and can be absorbed in the definition of 



66 



the vertex operators, and we find the scattering amplitudes at temperature T. For 
instance, 



S(Ei, E% E3) 



9str^{E\ + E2 — -E3) E1E2E3 1+ 



9str 



(10.20) 



24 



(1 + iE 3 )(2 + iE 3 ) (1 - i£ 3 + E\ + £ 2 2 + (2ttT) 2 ) + . . . 



Is the concept of finite-temperature S'-matrix well-defined? In the usual theories the 
answer is negative, because interaction with the heat bath stops a particle before it 
reaches the interaction region. However, in the hamiltonian (7.19) the interaction is 
exponentially localized in a finite region of space near r = 0. Thus, the interaction 
with the heat bath is negligible at infinity, and the concept of S'-matrix still appears 
to be meaningful. It would be interesting to compare the fermionic results like eq. 
(10.20) with bosonic finite temperature calculations with the hamiltonian (7.19). 

In conclusion, we should also mention the fascinating speculation due to Witten 
that the c = 1 matrix model describes physics in the background of a two-dimensional 
black hole. 51 * The black hole has finite temperature, and its Euclidean description 
is, perhaps, related to the compact c = 1 model with a definite radius R. It would 
be remarkable if the formulae reported in this section could be interpreted in terms 
of massless particle scattering off the black hole and in terms of Hawking radiation. 



★ For a detailed discussion, see H. Vcrlinde's lectures in this volume. 
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11. DISCRETIZED TARGET SPACE. 



In this section we will study another interesting model: the theory of surfaces 
embedded in a discretized real line with lattice spacing e. As shown in section 8, 
the lattice duality transformation relates it to the model of surfaces embedded in a 
circle of radius 1/e. Although for genus > this transformation generates 2h extra 
variables which violate the precise equivalence of the two theories, the basic physical 
effect - the Kosterlitz-Thouless phase transition - is of the same nature. 

In fact, the model with discretized target space is very interesting in its own right. 
Here we can test how string theory responds to introduction of a small lattice spacing 
in space-time. It has been argued 52 that in this respect string theory is very different 
from any field theory: if the lattice spacing is smaller than some critical value, then 
the theory on a lattice is precisely the same as the theory in the continuum. This is to 
be contrasted with any known field theory, where the continuum behavior may only 
be recovered as the lattice spacing is sent to zero. The exactly soluble c = 1 string 
theory is an ideal ground for testing this remarkable stringy effect. Below we give the 
exact solution of the model with the discretized target space and confirm that, for 
e < e c , the lattice is smoothed out so that the theory is identical to the embedding 
in R 1 . 

The matrix model representation of the partition function is now in terms of an 
integral over a chain of M matrices with nearest neighbor couplings 53 

M 



Z ^=flJ DN2 ®i™P -/3^^Tr($ i+1 -^) 2 + eTrW(^ 



;n.i) 



It is easy to check that the perturbative expansion of this integral gives the statistical 
sum of the form 



h A i=l rii (ij) 

which is precisely what is needed to describe the embedding in a discretized real line 
with lattice spacing e. As in the matrix quantum mechanics (2.6), the integral (11.1) 
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can be expressed in terms of the eigenvalues of the matrices $j. The only modification 
here is that, instead of quantum mechanics of N identical non-interacting fermions, 
we now find quantum mechanics with a discrete time step e. Thus, 



lim — —4- = ^ In m , 



(11.3) 



where are the N largest eigenvalues of the transfer matrix 



Vifi(x)= / dyK(x,y)fi(y), 



K(x,y) 



JL 

2vre 



1/2 



(11.4) 



cxp 



(3 Ux-yf 



+ e(W{x) + W(y)) 



Although it is hard to solve this problem exactly, as usual things simplify in the con- 
tinuum limit where we expect to find universal behavior controlled by the quadratic 
maximum of W(x). Indeed, if we take 



W(x) = -\x 2 + 0(x 3 ) 



11.5) 



and rescale x^ffi = z, we find 



K(z,w) = exp 



y/2 



7T6 



(z — w) 



■2, + rl z2+ro2+0 (^> (/+ro3) 



(11.6) 



Thus, the terms beyond the quadratic order are suppressed by powers of (3 and are, 
therefore, irrelevant. The quadratic problem can be solved exactly through finding the 
equivalent quantum mechanical problem with Planck constant 1//3' and hamiltonian 
H(e) such that 

K(x, y ) = (x\e-^' H ^\y) . (11.7) 



Then, ^ = exp(—ef3'ei) where ej are the N lowest eigenvalues of H(e). For the 
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quadratic transfer matrix, 



K(x,y) 



JL 

2ne 



1/2 



cxp 



(3 f(x - yf 



- \e(x 2 + y 1 



(11.8) 



„2 2 2 

the hamiltonian is also quadratic H = ^ — Comparing 



{x\ e 



-e/3'H 



\y) 



m \ V 2 
mujp \ 



2n sin we 



cxp 



with eq. (11.8), we find 



muif3' 



[(x 2 + y 2 ) cot cue — 2xysvn 1 uJt) 



(11.9) 



cos ecu = 1 — |e 2 , 



sin(c<Je) 



11.10) 



Thus, we have reduced the problem of random surfaces embedded in discretized target 
space to quantum mechanics of iV fermions moving in an inverted harmonic oscillator 
potential, i.e. to random surfaces embedded in R l \ Eq. (11.10) dictates how the cur- 
vature of the potential and (3' depend on e. Recalling that a' = 1/uj 2 and g st = ^rj^, 
we find that the target space scale and the coupling constant of the equivalent con- 
tinuum string theory depend on e. Thus, changing the lattice spacing in target space 
simply changes the parameters of the equivalent string theory embedded in continu- 
ous space. This effect was originally demonstrated in higher-dimensional string theory 
using completely different techniques. 52 It appears to be a general property of string 
theory. 

Intuitively, we do not expect the equivalence with the continuum string theory 
to be there for arbitrarily large e. Indeed, eq. (11.10) shows that, if e exceeds e c = 2 
then the solution for u becomes complex, which is a sign of instability. Since large 
lattice spacing corresponds, in the dual language, to small radius, it is clear that the 
instability is with respect to condensation of vortices. Thus, for e > e c , we expect 
the vortices to eliminate the massless field on the world sheet corresponding to the 
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embedding dimension, so that the pure gravity results. This is easy to check for very 
large e where the sites of the matrix chain become decoupled, 



lim 



lnZ(e) 
M 




D N 2 Q e -p e TrW(*)^ 



(11.11) 



and we obtain the one-matrix model well-known to describe the pure gravity. 

We will now give a continuum explanation of the Kosterlitz-Thouless transition 
from the small e phase, where the lattice is irrelevant, to the large e phase, where 
it has the dominant effect. We will replace the discretized real line by a continuous 
variable X with a periodic potential V(X) = J2 n> oa n cos(^-X), which implements 
the effects of the lattice. The conformal weight of the operator J d 2 a cos(^-X) is 

2 2 / 

h = - ^ Q . Thus, for small e all the operators perturbing the action have h > 1 
and are irrelevant. This is the essential reason why the string does not feel a lattice 
in target space with spacing < 0(Vc7). The phase transition to lattice-dominated 
phase takes place when the perturbation with n — 1 becomes relevant, i.e., for 



We can compare this result with the position of the K-T transition in the matrix 
chain model, where we found e c = 2 and y/cV(e c ) = 2/tt. The agreement of these 
values with the continuum relation (11.12) strengthens our explanation of the phase 
transition in the matrix model. It is satisfying that we have found the matrix model 
confirmation of a general effect, smoothing out of a target space lattice, which applies 
to string theories in any dimension. 
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12. CONCLUSION. 



I hope that I have convinced the reader that the 2-dimensional string theory is 
a highly non-trivial toy model for string theory in higher dimensions. Its matrix 
model formulation as a sum over surfaces embedded in 1 dimension is an example of 
a perfectly regularized generally covariant definition of the Polyakov path integral. 
It also turns out to be remarkably powerful, giving us the exact solution of a non- 
trivial string theory. Many of its physical features, such as the R — > 1/R duality 
(and its breaking due to the vortices), smoothing out of the target space discreteness, 
presence of poles in the correlation functions, etc., carry over to string theories in 
higher dimensions. Also, even though this string theory is two-dimensional, it has 
some interesting remnants of transverse excitations, manifested in the isolated states 
at integer momenta which generate the W^i+oo algebra. Thanks to the matrix models, 
we have acquired a wealth of exact information on the perturbative properties of 
this string theory. Unfortunately, we still do not have a detailed understanding or 
interpretation of most of it in the conventional continuum formalism. Hopefully, this 
gap will be closed in the not too distant future. 
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